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potential  distributions  are  taken  into  account.  Solutions  are  possible 
for  wide  ranges  of  the  values  of  tne  relevant  parameters,  and  are  applic- 
aole  to  the  sheaths  and  fluxes  of  emitting  spacecraft,  and  to  the  current- 
voltage  characteristics  of  emitting  Langmuir  probes.  The  formulation  for 
the  orbit  analysis  is  based  on  the  author's  (1973,  1975)  Turning-Point 
function  method,  as  opposed  to  the  metnod  of  the  Effective-Potential  func-  * 
tion  exemplified  by  tne  work  of  La  framboise  and  Bernstein  and  Rabinowitz. 

The  two  formulations  have  been  shown  by  the  autnor  to  be  analytically 
equivalent  and  to  produce  identical  numerical  results.  However,  in  the 
author's  opinion  the  Turning-Point  formulation  is  simpler  and  more  effi- 
cient than  the  Effective-Potential  formulation. 

Solutions  are  presented  for  a number  of  sample  problems,  using  a code 
PARKSS  based  on  the  analysis  of  this  report.  The  discussion  of  these 
solutions  includes  comparisons  with  previous  theories  and  experimental 
data.  In  these  problems  the  effects  are  considered  of  variations  of 
(a)  surface  potential,  (b)  monoenergeti c versus  Maxwellian  emission 
velocity  distributions,  (c)  emission  flux  relative  to  ambient  fluxes 
("'weak*  versus  "strong*  emission),  ana  (d)  emission  temperature 
relative  to  ambient  temperatures  (*colc*  versus  "warm''  emission),  on 
potential  barriers  and  net  surface  fluxes. 


1.  INTRODUCTION 


This  report  is  concerned  with  rigorous  calculations  or  the  structure 
of  the  charged-particle  sheath  about  a symmetric  body  immersed  in  a plasma 
and  emitting  electrons  from  its  surface.  The  body  can  be  a space  probe 
and  the  emission  Cu.n  consist  of  photoelectrons  or  secondary  electrons.  This 
becomes  an  important  effect  when  the  ambient  plasma  is  hot  and  tenuous,  say 
about  a spacecraft  outsice  the  plasmaspnere,  and  can  influence  electric 
field  measurements  (Grard  et  al , 1973)  or  low-energy  particle  experiments 
(Whipple  and  Parker,  1969).  The  model  to  be  discussed  is  equally  app 
ble  to  the  problem  of  an  electrostatic  Langmuir  probe  in  a laboratory  piisma 
where  the  emitted  electrons  are  thermionic  or  secondaries  (Chang  and 
BienKowski,  1970). 

Spherical ly-symmetric  sheaths  have  been  treated  rigorously  in  the 
absence  of  emission  by  Laframboise  (1966)  and  Parker  (1975).  The  effects 
of  emission  on  spherical  bodies  immersed  in  plasmas  have  been  treated 
approximately  by  a number  of  workers  (cf.  Whipple,  1976  and  references 
cited  therein),  but  very  little  rigorous  theoretical  work  has  Deen  reported. 
Schroder  (1973)  reports  tnat  ne  used  the  method  of  Laframboise  (1966)  to 
include  pnotoelectrons , but  gives  no  computational  details.  This  is  a sig- 
nificant omission  because  the  Laframboise  method  is  based  on  the  assumption 
tnat  tne  potential  profile  is  monotonic  whereas  Schroder's  results  involve 
nonmonotonic  potential  profiles  (potential  minima),  which  would  require  a 
significant  modification  such  as  that  of  the  present  report. 

There  is  ongoing  work  on  dynamic  time-simulation  calculations  for  tne 
spherical  sheath  problem  (P.  Rothwell , A.  Wilson,  private  communications, 
1976).  These  solutions  have  been  found  to  agree  (after  the  oscillations 
nave  become  sufficiently  small)  witn  the  author's  static  solutions  (Pur-er, 
1975  and  this  report)  at  suff 1 ciently  long  times.  However,  certain  ai**i- 
culties  are  encountered  with  dynamic  methods  which  are  absent  in  static 
methods.  For  example,  in  seeking  steady-state  solutions  one  nav  find  there 
are  two  or  more  different  time-scales  in  tne  problem,  which  can  make  time- 
simulation  relatively  expensive.  One  nay  also  inadvertently  introduce 


"numerical"  collisions  so  that  particles  oecome  inadvertently  trapped,  in 
a potential  minimum  for  example.  In  addition,  there  may  be  excessive  noise. 

A rigorous  spherically-symmetric  model  is  needed,  to  provide  relatively 
inexpensive  but  valuable  guidance  regarding  the  assessment  of  sheath  struc- 
tures in  3-dimensional  problems.  The  present  computer  model,  called  PARKSS 
(Parker  Spherical  Sheath),  has  been  developed  with  this  aim.  It  is  an  exten- 
sion of  the  author's  previously  developed  computer  model  (Parker,  1975),  but 
modified  to  take  any  nonmonotonic  potential  distribution  into  account.  It 
should  be  mentioned  that  a time-simulation  study  by  Soop  (1972)  for  a con- 
ducting sphere  illuminated  on  one  side  in  a vacuum  has  shown  that  t.i.  .sym- 
metry in  the  charge  distrioution  is  small.  Hence,  a spherically  symmetric 
model  should  be  a good  approximation  for  the  sheath  around  a spherical  space 
probe.  In  addition,  the  problem  of  a long  circular  cylinder  (which  has  clear 
implications  for  space,  e.g.,  instrument  booms,  as  well  as  laboratory  appli- 
cations) is  so  similar  to  that  of  a sphere  that  only  slight  modifications 
are  required.  Accordingly,  the  sphere  and  cylinder  are  treated  on  an  equal 
basis  in  the  analysis  of  this  report. 

In  Section  2 the  assumptions  and  basic  foundations  of  the  Vlasov  prob- 
lem are  discussed.  This  refers  to  the  evaluation  of  charged-particle  den- 
sities and  fluxes  with  arbitrary  potential  distributions.  In  Section  3,  the 
solution  of  Poisson's  equation,  and  the  Poisson-Vlasov  iteration  procedure 
to  obtain  mutual  consistency  of  the  distributions  of  potential  and  charged- 
particle  densities,  are  treated.  The  kinetic-theoretical  problems  of  the 
evaluation  of  densities  and  fluxes  are  discussed  in  Sections  4 and  5, 
respectively.  Here  we  consider  Maxwellian  and  monoenergetic  velocity  dis- 
tributions, for  ambient  and  emitted  particles,  and  for  spheres  and  cylinders. 
The  moment  integrals  are  represented  as  integrals  over  energy  and  angular 
momentum.  In  Section  6 we  consider  the  orbit  analysis  required  to  evaluate 
the  angular  momentum  integrals.  The  formulation  for  this  analysis  is  based 
on  the  Turning-Point  function  as  opposed  to  the  Effective  Potential.  The 
■Turning-Point  formulation  is  considered  to  be  much  simpler  in  practice  and 
more  efficient  than  tne  Effective-Potential  formulation.  The  equivalence 
between  the  two  formulations  has  been  proved  analytically  and  numerically 
(Parser,  1975),  but  is  here  demonstrated  in  one  of  several  examples  by 
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re-aeriving  on  tne  oasis  of  tne  Turning-Point  formulation,  the  formulas  derived 
earlier  by  Bernstein  and  Rabinowitz  (1959)  wno  used  the  Effective  Potential. 

Section  7 and  Apperdix  A are  concerned  witn  the  energy-quadrature 
method  wherein  the  moment  integrals  over  energy  are  approximated  by  finite 
sums  over  monoenergetic  factors.  The  method  of  Parker  (1975)  is  nere  moai- 
fied  to  deal  witn  nonmonotonic  potential  distributions  anc  arbitrary  energy 
ranges  wnich  arise  in  tne  emission  problem. 

In  Section  8,  a number  of  sample  spherical  solutions  obtaineo  with  the 
code  PARKSS  are  discussed.  The  effects  of  (a)  emitted  electrons  without 
the  ambient  plasma,  (b)  ambient  plasma  without  emission,  and  (c)  the  > 

ation  of  emitted  electrons  with  the  ambient  plasma,  are  compared.  When  there 
is  emission  there  is  generally  a negative  potential  minimum  (or  barrier) 
somewnere  Detween  the  probe  ano  infinity,  .ts  qualitative  behavior  is  as 
follows.  For  large  positive  surface  potentials,  the  minimum  has  negligible 
depth  and  occurs  at  a large  radius.  As  the  surface  potential  is  reducec, 
the  minimum  moves  inward  and  becomes  significantly  deeper.  At  a negative 
va'je  of  surface  potential,  the  minimum  occurs  at  tne  surface,  and  at  tne 
surface  the  field  is  therefore  zero.  With  more  negative  surface  potentials 
the  potential  profile  rises  monotonically  to  zero.  For  positive  surface 
potentials,  the  solutions  are  similar  to  tnose  presented  by  Schrbder  (1973). 
However,  tnere  are  some  differences  with  Schroder's  results.  It  is  shown 
also  that  emission  raises  the  equilibrium  potential,  as  expected.  The 
effects  of  strong  versus  weax  emission  (based  on  relative  fluxes),  and 
monoenergetic  versus  Maxwellian  distrioutions  of  emission  velocity,  are 
compared.  For  weax  and  cold  emission  from  a negative  oroDe,  the  emitted 
electrons  “run  downhill,"  ano  their  density  falls  rapidly  and  becomes  asymp- 
totically proportional  to  tne  inverse  square  of  the  radial  distance.  In  a 
small -Debye-number  probe  problem,  the  point  of  quasi  neutrality,  beyona  wnich 
positive  and  negative  charges  balance  to  within  a given  small  fraction,  is 
unchanged  oy  adding  electron  emission.  Rigorous  solutions  are  also  obtained 
for  a problem  treated  approximately  py  Whipple  (1976)  in  attempting  to  inter- 
pret ATS-6  spacecraft  data.  Tne  spherically-symmetric  model  was  aGoptec  to 
see  whetner  potential  barriers  indicateG  by  the  Gata  could  be  reproduced. 
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At  various  points  throughout  Section  8 equilibrium  potentials  are  indicated. 
The  limit  of  cold  emission  from  a negative  probe  is  derived  in  Appendix  B. 


The  computer  program  PARKSS  implementing  the  analysis  of  this  report 
is  presented  in  Appendix  C. 


2.  FOUNDATIONS 


The  spacecraft  is  modelled  by  a conducting  sphere  or  cylinder  while 
tiie  ambient  plasma  is  isotropic,  collisionless,  ana  stationary  with  respect 
to  the  spacecraft.  Pnotoelectrons  and  secondary  electrons  are  assumed  to 
be  emitted  isotropically  from  the  spacecraft  surface.  In  the  case  of  a 
Langmuir  probe  the  emission  can  consist  of  either  thermionic  electrons  or 
secondaries  due  to  bombardment.  Magnetic  field  effects  are  nenlected. 

Far  from  tne  probe  or  satellite,  the  distribution  functions  for  electrons 
and  ions  in  the  ambient  plasma  are  treated  in  tne  present  report  as  e. 
Maxwellians  or  monoer.ergetic.  The  same  applies  to  tne  distribution  func- 
tion for  emitted  electrons  at  the  surface.  The  problem  is  assumed  to  be 
independent  of  time.  It  is  desired  to  calculate  the  SDotial  distributions 
of  electric  potential  and  of  electron  and  ion  densities,  as  well  as  the 
currents  of  ions  and  electrons  exchanged  between,  tne  probe  anc  tne  envi- 
ronment. Tne  basic  problem  assumes  a fixed  potential  on  the  surface,  but 
it  is  a natural  extension  to  find  tne  floating  ootential  by  solving  a 
sequence  of  problems  until  current  balance  has  been  achieved. 

In  accord  with  the  foregoing  assumptions,  a particle  distribution 
function  f (which  is  the  density  in  phase  space)  satisfies  the  time- 
independent  Vlasov  equation, 


This  states  that  f,  which  is  a function  of  tne  vector  positior  R ana  or 
the  vector  velocity  V,  is  constant  along  an  orbit  in  an  electric  field  *' 
potential  e(R),  where  q and  m denote  respectively  tne  eta roe  ,-.rd  mass  a* 
tne  particle.  In  Eg.  ( 1 ) , . and  denote  gradient  o t 'ators  with  res  ect 
to  position  and  velocity  space,  respectively. 
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Assuming  one  lias  solved  L<j . (1)  for  f subject  to  the  appropriate 
boundary  conditions  (the  "Vlasov  problem"  - sec  Secs.  4-6),  one  may  then 

v w 

obtain  the  particle  density  f!(R)  at  a point  R by  integrating  f over  veloc- 
ity space: 

N(R)  - | ;[  f (R ,V)  dpif  (2) 

J J J 

where  dpV  represents  a p-dimensional  volume  element  of  the  velocity  space. 

For  the  spnerical  problem  p=3,  while  p-2  for  the  cylindrical  problem.  This 
integral  represents  the  zero-tn  moment  of  the  distribution.  Simile..'.  , one 
may  obtain  the  flux  or  current  density  at  a point  R by  defining  V as  the 
component  of  velocity  of  interest  and  integrating  the  product  fV'n  over 
velocity  space  (the  first  moment): 

j(R)  = jj]  f (R.V)  Vndpv  (3) 

Or.  a spnerical  or  cylindrical  surface,  V^  becomes  V , the  raaial  component 
of  velocity. 

One  may  determine  the  local  value  of  f(R,V)  at  an  arbitrary  point  R, 

—*•  ~0. 

by  considering  the  orbit  of  a particle  arriving  at  R with  the  velocity  V. 

If  the  orbit  is  traced  backwards  and  a point  is  eventually  reached  where 
f is  known  (namely,  the  "source"  at  the  surface  or  at  infinity),  then  tr.e 
local  value  of  f is  identical  to  tne  known  value  at  the  source.  It  the 
orbit  is  found  to  be  closed  upon  itself,  its  population  can  consist  on T. 
of  trapped  particles  which  can  nave  arrived  there  through  a coliisionai 
mechanism.  It  is  assumed  here  that  the  closed  orbits  if  any  are  unpopu- 
lated (that  is,  f=G  for  these).  This  assumption  is  usually  made  for 
tractability  purposes,  based  on  a plausibility  argument  vLafr.;. ' ooise , 

1966;  Whipple,  1976),  but  its  ultimate  justification  is  difficult  and 
requires  a quantitative  coliisionai  theory.  Such  a theory  is  available 
(Parxer , 1973) , but  nas  not  as  yet  been  applied  to  the  present  sheath  prob- 
lem. 

Since  we  have  assumed  isotropic  sources  at  the  surface  and  infirnt., 
ar.a  since  total  energy  is  conserved  along  an  orbit,  tne  velocity  distripjtions 
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can  dtpend  only  on  the  particle  energy  at  the  source.  Denoting  by  V the 
particle  speed  at  infinity,  and  by  Vs  the  particle  speed  at  the  surface,  we 
may  define  distinct  distribution  functions  for  ambient  and  emitted  pe- 
tioles, depending  on  the  distributions  at  the  respective  sources: 


l".Vj  F„(MV‘/2) 

(4 ) 

<-.»J  FJmvf/2) 

(5) 

<V5s>  - Fs(n,Vs2/2) 

(6) 

Here  the  subscripts  i_  and  e refer  to  the  ambient  ions  (of  mass  M)  and  elec- 
trons (of  mass  m),  respectively,  while  subscript  s refers  to  tne  surface- 
emitted  electrons.  The  arguments  of  F^  and  F$  are  relatea  to  the  total 
energy  which  is  constant  along  orbits.  Then  the  distribution  function  for 
the  ions  at  any  point  R can  be  written: 

f1  (R.V)  = 5i (R,V)  F„(Mv£/2)  v 7 ) 


For  tne  electrons,  there  are  separate  distribution  functions,  one  desc-io- 
ing  the  electrons  from  infinity  (F  ),  and  tne  other  those  from  tne  surface 


*e 


(R.V) 


fia(R,V)  F (mV“/2) 


£ 

I 


s 


(R.V) 


$S(R,V)  F$( 


o 


with  the  subscripts  e ana  s referring,  respectively,  to  tee  electrons 
infinity  (ambient),  ana  tnose  from  tne  surface.  Tne  "del tu-feccors " , 
(R,V)  contain  tne  orbit  information,  that  is,  regarnin.  wnether  tne 
orbit  connects  with  the  appropriate  source.  Thus,  a.  . : .a 

from  infinity,  while  5 ~0  if  the  ion  comes  fro;  the  p^obe  -rfae 


' , b 
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<5e=l  and  6^=0  if  the  electron  comes  from  infinity,  while  6^0  and  6S  = 1 if 
the  electron  comes  from  the  surface.  That  is,  a non-zero  value  for  6 sig- 
nifies that  the  orbit  is  "occupied."  Since  the  6 ' s are  step-functions,  and 
the  F-functions  are  known  (to  be  defined  in  Sec.  4),  it  is  only  necessary 
to  locate  the  boundaries  of  tne  regions  in  velocity  space  associated  with 
the  two  different  types  of  particle  sources.  For  tne  symmetric  problems 
of  interest  (spheres  and  cylinders),  the  regions  in  velocity  space  will  oe 
found  to  be  further  subdivided  according  to  types  of  orbits,  such  as  tnose 
which  contribute  once  due  to  a single  puss  tnrough  the  radius  of  i r.f-rest, 
and  those  which  contribute  twice  because  the  particle  reverses  direct-. u 
after  passing  the  radius  of  interest,  and  then  re-crosses  the  radius  of 
interest.  Thus,  the  6-factors  which  ao  not  vanish  can  take  on  the  value  2 
as  well  as  unity,  as  will  oe  shown  in  Sec.  6. 

Finally,  the  relation  between  tne  local  speed  V and  tne  speed  at  tne 
source  (V^  for  ambient  particles,  and  V$  for  emitted  particles)  is  given 
by  the  total-energy  relations 

Ambient  Ions 

KV» . 

-o-  -*•  eo  = -o-  + es  JO) 


Ambient  Electrons 


~2~  + ec>  = -jj—  + eoo 


Emitted  Electrons 


mV2  mVs\ 

__  + eo  = — + e$s 


wnere  « is  the  electric  potential,  cr  is  the  electric  potential  at  the  sur- 
face, and  is  the  potential  at  infinity  assumed  to  be  zero;  Eq.  (10/  -.s 
written  for  singly-charged  ions,  and  e is  the  signed  electron  cnarge. 
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3.  POISSON'S  EQUATION  (CALCULATION  OF  POTENTIALS) 

AND  ITERATION  PROCEDURE  FOR  POISSON- VLASOV  SOLUTION 


The  electric  potential  *5  is  obtained  from  the  solution  of  Poisson's 
equation,  which  may  be  expressed  in  terms  of  the  ambient-ion  number  density 

N . , one  ambient-electron  number  density  N , and  the  surface-electron  nu.oer 

1 e 

density  N : 

J s 

v2*  = 47re(Ns  + Ne  - \, ) (13) 

where  e denotes  the  rr.agnituce  of  tne  electron  charge.  It  is  convenient  to 
write  this  equation  in  dimensionless  form,  ana  in  terms  of  the  radial  coor- 
dinate appropriate  to  a sphere  or  cylinder,  where  the  radius  of  tne  sphere 
or  cylinder  is  denoted  Dy  r . In  non-di mens lunaii zing,  one  may  choose  tne 
ambient  parameters  fL  and  T as  basic,  where  N.  denotes  tne  ambient-plasma 
number  density  and  T Genotes  the  ambient  electron  temperature.  One  also 

C 

defines  N as  the  Surface-electron  number  density  evaluated  at  the  surface, 
so 

Tnen  for  the  sphere  Poisson's  equation  may  be  expressed  in  terms  of  the  fol- 
lowing dimensionless  variables: 

r = radial  coordinate  divided  by  r. 

5 = electric  potential  divided  by  kT^/e 

n.  = ambient-ion  density  divided  oy  N 

n = ambient-electron  density  divided  by  N 
e J 0 

nc  = surface-electron  density  divided  by  N 

■')  \i 

Aq  = ambient-electron  Debye  length  ( kT  /4rrN  eu)  ’■ 
divided  oy  r (that  is,  tne  Debye  number) 

The  resulting  equation  for  tne  spherical  problem  may  be  written: 


d2(r^)  _ 
~d?  " “ 


r . so 

-v  n fj — + n„ 

A^  L S N0  e 
D 


r 


The  form  of  the  left-hand  side  implies  that  this  equation  is  to  be  solved 
for  r$,  from  whicn  0 is  obtained  by  division  by  r.  Alternatively,  one  may 

non-dimensionalize  by  choosing  as  basic  the  emission  parameters  M ana  T , 

s ^ -* 

where  T$  is  the  effective  temperature  of  an  emitted  Maxwellian  distribution. 
Then  the  equation  may  be  expressed  in  terms  of  the  same  dimensionless  vari- 
ables except  the  following: 

<t>  = electric  potential  divided  by  kT^/e 

/ A ■ 

= surface-electron  Deoye  length  (kTs/4nNsoe  ) f- 
di video  by  rQ  (that  is,  the  Debye  number 
based  on  the  surface  emission) 


Then  the  equation  may  be  written: 


d2  ( r<j> ) 
dr2 


DS 


N N 

o o 

I n + j-  n - rj— 

[ s Nso  e "so 


I ■»  i 
< » « 


Equations  (14)  and  (15)  refer  to  tne  spherical  problem.  For  the  correspond- 
ing cylindrical  problem,  one  may  write  instead  of  Eq.  (14): 


s N 


(16) 


where  the  new  variable  u is  defined  by  u=ln  r.  Equation  (15),  oaseo  on  thr 
surface  emission  parameters , may  be  similarly  replaced. 

Equations  (14)-(16)  may  be  solved  by  differencing  them  on  a set  of  gria 
points,  in  r for  tne  sphere  and  in  u for  the  cylinder,  and  solving  tne  c ■ refer- 
ence equations.  The  reason  for  choosing  the  forms  of  Eqs.  (14)-(16)  is  that 
a simple  second-difference  operation  may  be  employed,  namely, 

( r^ ) -j _i  " 2(r'?)1  + (r<1),+-|  = (Ar)2(RHS) . i \ 

centered  at  the  i-tr.  grid  point  r.  for  the  spherical  problem,  and 

1 
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jt 

vlj_-j  _ + ^j+1  ~ (Au)  (RHS)j  ( IS) 

centered  at  the  j-th  grid  point  u^  for  the  cylindrical  problem.  In  Eqs. 

(17)  and  (18),  the  right-hand-side  factors  (RHS),-  and  (RHS).  denote,  re- 

* J 

spectively,  the  values  of  the  right-hand  sides  of  Eq.  (14)  and  Eq.  (16), 
for  example,  evaluated  at  the  center  grid  point;  and  Ar  and  au  denote  the 
interval  of  a uniform  grid. 

The  "double-sweep"  method  of  solution  for  such  tri -diagonal  systems  or 
linear  equations  is  well  Known  (Varga,  1962,  for  example).  Floating  Dcun- 
dary  conditions  to  represent  the  boundary  condition  at  infinity,  appro..  ■ .u 
for  electrostatic  probe  and  satellite  problems,  and  the  solution  for  non- 
uniform  grids,  are  treated  in  Parker  (1975). 

Given  the  method  for  evaluating  the  ion  and  electron  densities  (to  oe 
discussed  in  the  next  section),  the  numerical  solution  of  the  sheath  prob- 
lem requires  that  of  two  fundamental  sub-problems,  the  "Poisson  ProblenV 
(this  section)  and  the  "Vlasov  Problem"  (next  three  sections).  The  Poisson 
Proolem  is  solved  to  yield  the  potentials  at  the  grid  points  when  tne 
charged-particle  densities  are  given,  and  conversely,  the  Vlasov  Problem  is 
solved  to  give  tne  charged-particle  densities  at  the  grid  points  when  the 
potentials  are  given.  For  mutual  consistency  an  iterative  process  is  re- 
quired. 

A suitable  procedure  for  tne  iteration  is  given  by  Parker  and  Sullivan 
(1974)  and  Parker  (1975).  In  tne  procedure,  one  mixes  successive  iterates  or 
tne  charge  density  (or  alternatively  the  potential  iterates).  This  means  that 
one  mixes  a fraction  (a)  of  the  most  recently  computed  charge-density  distri- 
bution (RHS  in  Eqs.  (17)  and  (18))  with  the  complementary  fraction  (1-0  of 
the  previously  used  charge-density  distribution.  For  small  Debye  numbers, 
the  fraction  a must  be  small.  There  is  an  upper  limit  for  «,  and  correspond- 
ingly a minimum  number  of  iterations,  such  that  convergence  is  assured. 

These  depend  on  the  type  of  boundary  condition  used  to  represent  infinity, 
as  well  as  on  the  number  of  Debye  lengths  between  tne  surface  and  tne  outer 
grid  boundary.  It  can  be  shown  (a)  that  for  small  Deoye  numbers  th<  nmber 
of  iterations  required  is  proportional  to  tne  square  of  tne  inverse  Deb/e 
number,  and  (b)  that  no  mixing  is  required  for  convergence  if  the  boundary 
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of  the  grid  is  less  than  about  it  Debye  lengths  from  the  surface  (Par«er  and 
Sullivan,  1974). 


The  boundary  condition  where  the  potential  floats  in  accord  with  an 
assumed  inverse-square-law  behavior  is  physically  appropriate  (Laframboise, 
1966)  and  is  found  to  be  computationally  efficient  for  the  purely-am.Dient- 
plasma  problem  (Parker  and  Sullivan,  1974).  This  condition  is  also  oeir.g 
used  provisionally  in  the  present  program. 


L. 


4.  THE  VLASOV  PROBLEM:  CHARGED-PARTiClE  DENSITIES 


We  first  consider  Maxwellian  distributions , of  emitted  and  ambient 
particles,  for  the  spherical  problem  (subsections  A1  ana  A?)  and  for  the 
cylindrical  problem  (subsections  81  and  52).  Then  we  consider  correspond- 
ing monoenergetic  distributions  (subsections  C and  D).  Ambient-particle 
examples  have  beer,  previously  derived  by  Laframboise  (1966)  and  Parker 
(1973  and  1975). 

A.  MAXWELLIAN  DISTRIBUTIONS- SPHERE 

We  will  derive  in  detail  tne  number-density  integral  in  the  spherical 
problem  for  emitted  particles.  The  ambient-particle  number  density  formula 
will  be  shown  to  fellow  easily  from  tnat  of  the  emitted  particles. 

Ai . Maxwellian  Emitted  Particles,  Spnericu'  Problem 

The  number-density  integral  for  eni tteo  particles  at  a radial  position 
R in  the  spheri cal  problem  is  conveniently  written  as  a specialization  of 
tne  velocity  space  in  Eq.  (2)  with  p=3  to  spherical-polar  velocity  coordin- 
ates V,  6,  and  3,  where  V is  the  local  speed,  e is  the  local  polar  angle 
with  respect  to  the  radial  direction,  and  3 is  the  local  azimuthal  angle. 

3^  2 

The  velocity-space  volume  element  is  d V = V dV  d(cose)d3,  but  because  of 
tne  assumed  isotropy  the  integration  over  3 is  trivial  (yielding  the  factor 
2tt ) , and  Eq.  (2)  in  view  of  Eqs.  (4)  and  (7)  may  be  written: 

NS(R)  --  2 it  j F$(mV2/2)V2dV  j _ d-c(cose)  (19) 

(sphere)  "u 

wnere  <5  denotes  tne  "delta-factor"  for  tne  emitted  particles  (-6  in  Sec.  2 
but  here  tne  subscript  is  dropped),  which  depends  on  tne  nature  of  the 
There  are  in  general  three  distinct  ranges  of  integration.  Tnose  are 
a teci  with  tnree  correspondin'.:  types  of  contributing  orbits,  "Type  ' , " 

2,"  and  "Type  3"  as  are  defined  in  Sec.  6.  We  defer  further  discus; in  of 
tne  orbits  until  Sec.  6. 

For  a Maxwellian  distribution  of  emitted  particles  at  the  _.ur f.ic 

have 
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(20) 


4 


Umvf/2)  = ( 2 ) N 


3/2 


so(2ukT  j 


exp(-rnV^/2kTs) 


where  T&  is  the  Maxwellian  emission  temperature.  The  first  factor  of  2 in 
parentheses  represents  the  doubling  appropriate  to  a half-Maxwel lian  dis- 
tribution at  the  surface.  This  is  consistent  with  Chang  ana  Bienkowski 
(1970)  and  Schroder  (1973).  At  zero  field,  the  angular  integral  in  Eg. 

(19)  has  contributions  from  only  a hemisphere  of  directions  at  the  surface 
(involving  outgoing  particles  only).  This  situation  is  in  contrast  with  the 
ambient  plasma  in  which  there  is  a full  Maxwellian  at  infinity  (involving 
both  outgoing  and  incoming  particles,  at  zero  field);  in  the  ambient  dis- 
tribution corresponding  to  Eq.  (20)  the  factor  2 is  replaced  by  unity. 


Alternatively,  one  might  wish  to  consider  the  surface  distribution  to 
be  a full  Maxwellian,  but  in  the  absence  of  an  infinite  external  potential 
barrier  there  would  be  an  imbalance  of  outgoing  and  incoming  particles  just 
outside  the  surface,  resulting  in  a discontinuous  density  function.  This 
situation  could  probably  be  made  consistent  by  applying  suitable  re-definitions, 
but  such  a procedure  seems  neither  immediately  obvious  nor  physically  satis- 
fying. In  any  case,  the  normalization  chosen  affects  tne  relationship  between 
the  current  and  the  density.  Tne  factor  2 in  parentheses  will  be  retained 
here  for  tne  emitted  particles  but  will  be  deleted  for  the  ambient  particles. 

With  Eq.  (20)  for  F , Eq.  (19)  becomes 


MS(R) 


(2)211 


so 


(JL_ 

'2kT 


3/2 


Jo 


exp(-mV^/2kTs)V2dV 


6'd(cose ) (21) 


where  V$  is  connected  to  V by  Eq.  (12).  It  is  desirable  next  to  make  use  of 
the  symmetry  of  the  problem  to  convert  the  variables  of  integration  V and  u 
to  new  variables,  namely,  total  energy  and  angular  momentum,  botn  of  wnich 
are  conserved  along  the  orbits.  It  is  also  desirable  to  nondimensionalize 
tne  variables.  We  write: 
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r 

II 

50 

-S 

o 

n 

s 

= VMso 

0 

v" 

= mV2/2kT  , 

V d 

= mV2/ 2k' 

s 

s 

s 

9 

= es/xT^, 

^s 

= e*t/kT, 

(22) 


and  convert  to  dimensionless  total  energy,  namely,  Eq.  (12)  nondimension- 
alized  through  division  oy  kf  : 


2 2 
E =v  + $ = vs  + 


(23) 


and  the  square  of  the  angular  momentum,  nondimensionalizea  through  division 
by  2mr02kT  : 


,2  _ R2  V2 

0 " o 1 


s r'/si'r'6 


7 C 


For  the  sphere  the  angular  momentum  is  that  about  the  center  of  the  snnere, 
while  for  the  cylinder  (treated  next)  it  is  that  about  the  axis.  Tie  advar,- 
tage  of  converting  to  the  new  velocity-space  variables  (E,J  ) is  that  botr, 
are  constants  of  the  motion  and  completely  characterize  an  orbit  independent 
of  tne  variation  of  r along  that  orbit,  in  an  orthogonal  r-E-J  space  the 
oroits  are  thus  straight  lines  parallel  to  the  r-axis  (Parker,  1975),  anc 
the  orbit  analysis  is  simpiifiea. 

O 

Using  Eqs.  (23)  and  (24)  to  eliminate  v and  e in  favor  of  E anc  3^,  we 
obtain  from  Eg.  (21 ) : 


n (r)  = __  exp(y  ) exp(-E)dt •M  ( E) 
s /—  s ns 

/ . /Tt  Jo 

(sphere) 

where  T1 (E)  is  tne  spherical  morioenergeti c contribution  to  r,  a^e  to  t-r.rvv 
ns 

E.  This  monoenergetic  factor  is  given  by 
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( 5 to  be  evaluated  in  Sec.  6) 


2 

which  is  comprised  of  the  integral  over  J with  E fixed. 

A2.  Maxwellian  Ambient  Particles,  Spherical  Problem 

The  Maxwellian  number-density  integral  for  ambient  particles  at  a 
dimensionless  radial  position  r in  the  spheri cal  problem  is  so  similar  to 
that  of  the  emitted  particles  that  the  final  expression  can  be  obtained 
readily  from  Eqs.  (25)  and  (26).  Thus,  we  set  ^c=0  since  it  now  represents 
the  potential  at  infinity,  delete  the  factor  2 in  parentheses,  and  cotain, 
for  either  the  ambient  ions  or  the  ambient  electrons: 

nanb(r>  = 2 ["  exp(-E)dE-Hns(E)  (27) 

(sphere)  7r  ^o 

where  M (E)  is  the  monoenergetic  contribution  of  the  energy  E.  This  fac- 
tor is 

I rv  E-<j> 

Mn$(E)  = £ ; 6*d[E-<j>-J  /r  ] (23) 

(sphere)  o 

(s  to  be  evaluated  in  Sec.  6) 
where  we  have  defined: 

namb  = ‘^i  e^o  <:°r  ions’  electrons 
<?  = eo/kT.  for  ions,  electrons 

1 9 L 

T.j.T  = arabient-ion,  ambient-electron,  temperature 

E = (MV2/2  + |ej*)/kT.  for  ions  v2'5) 

1 (mV‘/2  - je|4>)/kT(i  for  electrons 

(K^/r^) (MV^sin^o )/2kT  for  ions 
o i 

^ - (A  /r^  / (lii'i/^ sir^O )/2kTL,  for  electrons 


j 


MAXWELL  LAN  DISTRI  BIT  IOMS-CYl  I NDER 


j. 

We  next  consider  Maxwellian  distributions , of  emitted  and  ambient  par- 
ticles, in  tne  cylindrical  problem. 

B1 . Maxwell; an  Emitted  Particles,  Cylindrical  Problem 

The  number-density  integral  for  eral tted  particles  at  a radial  position 

R in  the  cyl 1 ndrical  problem  is  obtained  when  one  specializes  the  velocity 

space  in  Eq.  (2)  witn  p-2  to  cylindrical  velocity  coordinates  V ana  6-  Here 

V is  the  speed  in  the  transverse  plane  perpendicular  to  the  cylinder  ..w 

and  0 is  the  (azimuthal)  angle  in  the  transverse  plane,  with  respect  to  tr, v_ 

2' 

radial  direction.  The  velocity-space  volume  element  is  d V=VdVde,  and 
Eq.  (2)  now  becomes 

N.vR)  = 2 ; Fe(mV?/2kTjVdV  " <5*cei  ^30) 

->  j b b b j 

(cylinder)  0 o 

v,.,ere  the  "delta-factor"  is  identical  to  that  of  Eq.  (19)  for  tne  sphere. 

For  a Maxwellian  distribution  of  emitted  particles  at  the  surface  we  nave: 

Fs(mV2s/2)  = (2)  Nso(^r ) exp(-mV2/2kTs)  (31) 

Tne  considerations  of  a half  Maxwellian  versus  a full  Maxwellian  discussed 
following  Eq.  (20)  apply  equally  well  here.  Tne  factor  2 in  parentheses 
will  be  retained  here  for  the  emitted  particles  rut  will  be  deleted  for  tne 
ambient  particles.  Equation  (30),  with  Eq.  (31),  may  do  converted  to  2-." 
space  witn  variables  nondimens i onal i zeo  as  in  Eqs.  (22)-( 24),  ana  we  .-rtair. 

n$(r)  = exp(qs)  ; exp(-E)dE*Mnc(E)  (3_ 

(cylinder)  c 

wr‘ere  Mnc^  is  tno  cy^ i ndrical  monocnergetic  contribution  to  n due  to 
enorqy  E. 

This  factor  is  given  by 


(33) 


V£) 

(cylinder) 


(2) 

If  J 
J0 


( 2 ) rT'  - 1 r 2 2.  , 

= — . 6 *d  (sin  lJ  /r"vE-$)]  > 


(s  tc  be  evaluated  in  Sec.  6) 


which  is  comprised  of  the  integral  oyer  o'"  with  E fixed. 


B2 . Maxwellian  Ambient  Particles,  Cylindrical  Problem 

The  Maxwellian  number-density  integral  for  ambient  particles  at  r in 
the  cvl i nqrlcai  problem  is  obtained  readily  from  Eqs.  (32)  and  (33).  Again, 
setting  <?s=0  ana  deleting  the  factor  2 in  parentheses  as  in  the  spnerical 
problem,  we  obtain,  for  either  ambient  ions  or  ambient  electrons: 


namh(r)  = j exp(-E)dE-Mnc(E)  (34) 

(cylinder)  0 

where  Mnc(E)  is  the  monoenergetic  contribution  given  by 

Mnc(E)  = 7 j 6*d{sin" 1 CO^/r1- ( E-p ) ] } (35) 

(cylinder)  0 

The  dimensionless  variables  are  defined  by  £q.  (29). 


C.  MONO ENERGETIC  QISTRIBUTIONS-SPHERE 

It  is  sometimes  convenient  to  approximate  a Maxwellian  distribution 
oy  an  equivalent  monoenergetic  one  (Laframooise,  1966)  and  for  this  reason 
the  monoenergetic  formulas  are  of  interest. 

Cl.  Monoenergetic  Emitted  Particles,  Spherical  Problem 

For  emi tted  particles  in  the  spheri cal  problem,  it  can  be  shown  tnat 
the  appropriate  monoenergetic  distribution  is 


* 


Fs(mv‘/2)  = 


(2)N 


3/1 


0O-mVj/2Ko) 


(36) 


where  D(x)  is  the  Dirac  delta-function,  and  K is  the  singular  value  of  the 

0 

kinetic  energy  at  the  source.  Substituting  Eq.  (36),  in  place  of  the 
Maxwellian  Eq.  (20),  into  Eq.  (19),  and /londimensional i zi ng  as  in  SuDsection 
A1  (Eqs.  ( 22 ) - ( 24 ) ) but  with  tne  normalizing  energy  kT$  replaced  by  K^, 
yields: 


n$  (sphere,  monoenergetic) 


. (2)  r/iTV* 

Jo 


5-d[ 1+v  -o-J2/r2j 


Lh 


(37) 


Inis  expression  becomes  identical  to  the  monoenergetic  factor  Eq.  (26)  when 
E is  replaced  by  l+$s. 

C2.  Monoenergetic  Ambient  Particles,  Spnerical  Problem 

For  ambient,  particles  in  the  spherical  proDlem,  tne  monoenergetic  num- 
ber density  is  obtained  readily  from  Eq.  (37),  oy  setting  ^=0  and  deleting 
the  factor  2 in  parentheses.  Thus: 

namb  (sphere,  monoenergetic) 

1 'fl-*  2 2l/2  . 

= o ■ 6‘d[l-v-u  /r  J ( 33 j 

Jo 

where  the  dimensionless  variables  are  defined  oy  Eq.  (29)  but  with  kT_  or 

kl  replaced  by  K . 
e r o 

D.  MONOENERGETIC  DISTRIBUTIONS-CYLINDZR 

In  contrast  to  the  spherical  problem,  it  is  doubtful  if  in  the  cylin- 
drical problem  it  makes  sense  to  approximate  a Maxwellian  distribution  :y 
tin  equivalent  monoenergetic  distribution  confined  to  the  tr,ir  verse  pi  >.nt -. 
..owevor,  the  appropriate  formulas  are  included  here  for  completeness,  with- 
out derivation. 


•i 


ns  ( cylinder,  monoenergetic) 


_ -I  2 Z /;> 

6*disin  [J  /r‘"(l+.?s-v>)  j ) 


n ^ ^linder’  monoenergetic) 


*d{  si  n~ 1 r^{ i-$)]  > 


where  all  quantities  have  been  previously  defined. 


THE  VLASOV  PROBLEM: 


5.  THE  VLASOV  PROBLEM:  CHAkGLO-PARTIClE  CuRRlNTS  (FluXES) 

The  oraer  in  whicn  tne  various  types  of  currents  (fluxes)  are  con- 
sidered will  follow  that  of  the  previous  section  on  number  densities. 
Namely,  we  treat  'axwellian  distributions  of  emitted  and  ambient  par- 
ticles for  the-  spherical  proolem,  ana  then  the  same  for  the  cylinorical 
problem.  Following  these,  the  corresponding  monoenergetic  distributions 
will  be  considered. 

A.  MAXWELLIAN  DISTRIBUTIONS-SPriERE 

We  derive  here  tne  flux  of  emitted  particles  in  the  soherical  problem. 
The  ambient-particle  flux  will  follow  readily. 

A'i . Maxwellian  Emitted  Particles,  Spherical  Problem 

The  current  or  flux  integral  for  emit tea  particles  in  tne  spnerical 
problem  is  given,  in  accord  with  Eq.  (3)  with  p=3,  by  multiplying  the  inte- 
g.  .nd  of  Eq.  (19)  by  Vr?Vcosa,  the  radial  component  of  velocity.  Thus, 
tne  emitted  flux  may  be  written: 


js (sphere)  = 


" , Fs(mVV2)VJoy  o-d(sin“o) 
J0  •‘o 


Here,  the  range  of  9-integration  is  reduced  by  one-half,  namely  (o,r/2), 
since  rlux  contributions  occur  only  over  a hemisphere  of  directions,  with 
given  by  Eq.  (20)  for  the  Maxwellian,  the  analysis  proceeds  in  a manner 
similar  to  that  for  the  density.  The  results  are  as  follows: 

jc(spnere,  Maxwellian) 

fwT 

= exP(0  exp(-E)dE*M  (E)  v 42 

i'  J ^ 


wnere  M,  (E)  is  tne  monoenergetic  contribution  given  by 

si  ^ 


(43) 


M.  (E)  = ^c-<J>^6*d(J2/r2) 

:s  > o 

(where  r-1  at  the  surface) 

where  all  dimensionless  quantities  are  defined  in  Sec.  4.  7he  neutral 
flux  (at  zero  field)  is  (2)Nso/kTs/3^m. 

A2.  Maxwellian  Ambient  Particles,  Spherical  Problem 

As  in  Sec.  4,  we  obtain  the  ambient-particle  expression  from  ::.e 
enitted-particle  expression  by  (a)  setting  $s=0,  (b)  deleting  tne  facto; 
in  Darentheses,  and  (c)  replacing  tne  emission  parameters  oy  the  corres- 
ponding ambient  parameters.  Thus  we  nave 

^amb^sphere?  MdXwellian) 

/kT . / 2 — M 
v 

= N / 

° i 

^ AT  e/2irm 

where  the  monoenergetic  factor  M-  (E)  is  identical  to  tnat  given  tr.  i...  (43), 
and  the  upoer  and  lower  square-root  factors  belong,  respecti vc’y , to  tr.o 
ambient  ions  and  electrons.  The  neutral  flux  (at  zero  field)  is  ’« . /kT/2-1  . 

3.  MAXWELL  I A, T DISTRIBUTIONS-CYLINDER 

31 . Maxwellian  Emitted  Particles,  Cylindrical  Proble- 

For  the  cylinder,  we  multiply  the  integrand  of  Eq.  (30)  by  Vv.  Vcor-, 
the  radial  comoonent  of  velocity.  Then  the  cylinder  emitted  flux  may  be 
wri tten: 


exp(-E)dE-M  (E) 
Jo 


f”  ? o n 

j (cylinder)  = 2 Fc  (niV‘72)VwdV  6*d(sine)  v45) 

s t s s 

;o  Jo 

where  the  range  of  o-integration  is  again  halved  to  reflect  tr.--  restriction 
ot  directions  to  tnose  contained  in  i nemicyi  moer.  With  F.  river,  oy  Eq. 

( 31 ) , we  obtai n: 


(cylinder,  Maxwellian) 


/ — S / 

= (2)N  / —5 — exp( p ) , exp(-E)d£-M.  (E) 

S0V  /m  s Jq 


wnere  the  monoenergetic  contribution  M-  (t)  is  given  by 

J 


M . ( E ) = o*d( ✓62/r2) 
Jc  Jo 

(where  r=l  at  the  surface) 


The  neutral  flux  (at  zero  field)  is  (2)hso/kTs/2nm. 


32.  Maxwell  1 an  Ambient  Particles,  Cylindrical  Problem 

We  obtain  the  ambient-oartic'ie  expression  from  the  emitted- particle 
expression  as  in  Sec.  A2  above.  This  yields: 


jamD (cylinder,  Maxwellian ) 


/2kT . / it  M 


!kV’ n J 


• : exp(-E)dE-n  (E) 

Jo  J 


where  we  take  the  upper  square  root  for  the  ions  and  the  lower  for  the  elec- 
trons; M.  (E)  is  given  by  Eq.  (47).  Tne  neutral  flux  (at  zero  field)  is 

h v^T/EVrrT. 

0 


C.  MONOENERGETIC  DISTRIBUTIONS- SPHERE 


Using  the  monoenergetic  oistribution  Eq.  (3o)  t,  Eq.  (41)  os. 
tne  following  monoenergetic  fluxes  for  tne  spherical  problem: 


j (sphere,  monoenergetic) 


= (2)N 


■0+<?s-<?) 


S*d(J2/r2 ) 


(where  r=l  at  the  surface) 


(49) 


for  the  emitted  particles,  ana  therefore 


^amb^spi',ere’  mon°energetic) 


( yftjm  ) 


..  j [ rO-9) 

= N0  { ) • ! 


d(u2/r2) 


/K0/8m 


(where  r=l  at  the  surface) 


(50) 


"or  the  ambient  ions  and  electrons,  respectively.  The  neutral  fluxes  (at 
zero  field)  are  ( 2 )NS0/KQ/ 8n;  and  NQ/Ko/8m,  respectively.  Note  that  for  the 
sphere  tne  monoenergetic  neutral  flux  is  equal  to  the  correspond! r.g  Maxwel- 
lian neutral  flux  if  Ko=4kT/-rr. 

D.  MONOENERGETIC  DIS7RI3UTI0NS-CYLINDER 

lie  reoeat  nere  tne  statement  of  Sec.  4D  that  it  seems  to  make  no  sense 
in  the  cylindrical  problem  to  use  a monoenergetic  distribution  in  the  trans- 
verse plane.  However,  the  formulas  are  included  for  completeness. 


j s ( cy 1 i nee r , monoenerge 1 1 c ) 


(where  r=l  at  the  surface) 


(51) 


"or  tne  e~i tted  particles. 


ana 


Jamb(cy‘inder*  rionoerier9et'lc) 

( / 2 'l 

j WJ™  : 

= N l ')•  *“’Vd(/J2/r2)  (52) 

0 ! / r : 

i /2V* ra  I 

(where  r=l  at  the  surface) 

for  the  ambient  ions  and  electrons.  For  the  cylinder,  the  monoenergeti c 
neutral  flux  becomes  equal  to  the  correspond! ng  Maxwellian  neutrai  flux  if 
K0=*kT/4. 


6.  THE  VLASOV  PROBLEM:  ORBIT  ANALYSIS 


In  Secs.  4 and  5 the  monoenergetic  contributions  involved  ranges  of 
integration  over  angular  momentum  represented  by  "delta-factors"  associated 
with  these  ranges.  The  delta-factors  will  be  defined  here  in  terms  of  the 
contriDutions  of  three  types  of  orbits,  namely.  Types  1,  2,  and  3 (Parker, 

1973,  1975),  as  illustrated  in  Fig.  1. 

Type  1 

The  Type-1  orbit  passes  from  infinity  to  tne  sphere,  one  way,  or  fro... 
tne  sphere  to  infinity,  one  way,  with  no  turning-point.  The  ingoing  orbit 
is  populated  by  a particle  from  the  ambient  plasma,  while  the  outgoing  orbit 
is  populated  by  a particle  emitted  from  tne  surface.  For  Tyre-1  orbits  tne 
delta-factor  is  unity. 

Type  2 

The  Type-2  orbit  starts  at  infinity  and  passes  inward,  to  a minimum 
raaius  (radial  turning-point)  outside  the  sphere.  It  subsequently  goes  out 
again.  These  are  two-way  orbits,  occupied  by  ambient  particles  going  in 
either  direction.  For  Type-2  orbits,  the  delta-factor  is  2.  (It  is  of 
course  zero  if  the  turning-point  is  outside  the  radius  of  interest.) 

IZSeJ. 

The  Type-3  orbit  starts  at  the  sphere  surface  and  passes  outward  to  a 
maximum  radius  (radial  turning-point),  subsequently  returning  to  the  spr.ere. 
These  orbits,  similar  and  inverse  to  the  Type-2  orbits,  are  also  two-way 
orbits;  they  are  occupied  by  emitted  particles  going  in  either  direction. 

For  Type-3  orbits  tne  delta-factor  is  2.  (It  is  of  course  zero  if  the  turning- 
point  is  inside  the  radius  of  interest.) 

Type  4 

This  closed  orbit  is  populatable  by  collisions.  Examples  are  giver  by 

Parker  (1973). 


1 


Consider  a typical  monoenergetic  integral  of  Secs.  4 or  5.  Denote 
its  value  by  M(E).  Tnen  it  nas  the  form: 

f o o 

M(E)=(const)  | <5  •dG(J‘'/r‘")  (53) 

J 

where  G(X)  is  a function  having  various  forms: 

G(X)=/K-X  spnerical  density 
G(X)=sin"  1 A/R  cylindrical  density 
G(X)=X  spherical  flux 
G(X)  = v/^’ cyl indrical  flux 

2 2 

where  the  total  range  of  X is  (0,K).  Here,  X represents  J /r  while  K 
represents  E-$. 

The  analysis  of  Parker  (1973)  snows  that  it  is  poss.ble  to  nave  Type-1 

and  Type-2  orbits  contributing  simultaneously,  or  it  is  possible  to  have 

Type-1  and  Type-3  orbits  contributing  simultaneously,  but  it  is  not  possible 

to  have  Type-2  and  Type-3  orbits  contributing  simultaneously.  Moreover,  if 

there  are  any  contributions  at  all,  at  least  some  of  these  come  from  the 
2 

vicinity  of  J =0  in  the  form  of  Type-1  orbits.  Hence,  if  the  point  r can  be 

2 

reached  energetically,  there  is  always  a lower  range  of  J for  whicn  tnere 

2 

are  only  Type-1  orbits,  and  there  may  be  an  upper  range  of  J for  which 
there  are  either  Type-2  or  Type-3  orbits.  Let  J-.2,  J02,  and  J2  denote 

2 ~ ??  ^ o o J 

critical  values  of  J as  follows  (with  0<J-|  or  0<J.j  <J^): 

2 

(0,J-  )=Kange  of  Type-1  orbits 

(J.j  ,j£  )=Range  of  Type-2  orbits  (54) 

2 2 v 

(J,  ,J,  )=Range  of  Type-3  orbits 
* 

At  the  surface,  the  flux  and  density  can  nave  contributions  only  fro;, 

Type-1  orbits. 

2 2 2 

In  terms  of  J^,  , and  J.  , t 

may  now  be  written  for  dens  it1': 


_ 7 0 


he  monoenergetic  integral  of  Eq.  (53) 


(55) 


1 


M(E)  = (const)  ;G(0,2/r2)-G(0)+2G(o//r^)-2G(J  2/r2), 


= (const)  i-G(0)+2G(J[?12/r‘1')-G(J1"/r2) 


wnere  m=2  (Type-2  oroits)  for  ambient  particles,  and  m=3  (Type-3  orbits) 
for  emitted  particles.  Tiie  aosolute-value  signs  take  into  account  tnat  G 
may  be  either  an  increasing  or  decreasing  function.  If  there  are  only 
Type-1  orbits,  tnen  Eg.  (55)  is  replaced  by 

M(E)  = (const) *G(J12/r2)-G(0); 


inis  equation  always  applies  to  fluxes.  It  may  be  obtained  from  Eq.  (55) 

2 2 

by  setting  J . 

J 3 m i 

2 

Next,  we  consider  how  the  critical  values  of  J are  determined. 

We  write  the  equation  for  the  conservation  of  energy  in  tne  dimension- 
less form 


E 


v 2 + 02/r2 
r 


(57) 


where  denotes  the  radial  component  of  velocity  (in  either  tne  spherical 

or  the  cylindrical  problem).  Then  along  an  orbit  with  fixed  total  energy  E 

2 

and  fixed  angu lar-momentum-squareh  J,  tne  radial  velocity  v^  will  remain 
finite  and  will  not  vanish  or  change  sign  as  long  as  the  following  inequal- 
ity is  satisfied: 


b = b 


+ J2/r2 


(53) 


Alternatively,  the  inequality  may  be  expressed  by: 

,2  2,r-  . 

J < g = r (E-<p) 


(59) 


Equation  (53)  may  be  analyzcG  to  find  tne  maxima  in  the  function  v,  which 
denotes  the  "effective  potential."  Alternatively,  Eq.  (59)  may  be  analyzed 
to  find  the  minima  in  the  function  g,  wmen  denotes  the  "turning-point 
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function."  Tne  analyses  of  Eqs.  (58)  and  (59)  to  classify  the  orbits  repre- 
sent, respectively,  tne  methods  of  the  Effective-Potential  Formulation  and 
the  Turni ng-Poi nt  Formulation.  The  effective-potential  approach  is  repre- 
sented by  the  works  of  Bernstein  and  Rabinowitz  (1959)  and  Lafrainboise  (1966), 
while  the  turning-point  approach  is  exemplified  by  the  works  of  Bohn  e al 
(1949),  Allen  et  al  (1957),  fledicus  (1961),  and  Parker  (1973,  1975).  Tne 
two  approaches  are  of  course  equivalent  (Parker,  1973,  1975).  They  involve 

different  projections  of  the  same  3-dimensional  phase  space  of  the  variables 
2 

r,  t,  and  J . From  the  definition  that  vr=o  defines  a turning-point,  Eq. 

(57)  with  vr=o  defines  a 3-dimensional  surface  representing  the  locu„  of 

turning-DOints.  A maximum  of  occurring  in  an  r-E  projection  (at  constant 

2 2 
J ) corresponds  to  a minimum  of  g occurring  in  an  r-J  projection  (at  constant 

E)  (Parker,  1975). 

The  Turning-Point  Formulation  is  considered  by  the  author  to  be  much 
simpler  in  practice  than  the  Effective-Potential  Formulation,  and  has  been 
shown  to  produce  identical  results  (Parker,  1975).  In  using  the  Turning- 
Point  Formulation  to  evaluate  monoenergeti c contributions  (that  is,  contri- 
butions to  the  density  or  flux  due  to  a given  energy),  one  explores  tne 
g-function  of  Eq.  (59)  for  least  values  (or  minima)  occurring  in  tne  radial 

range  r=l  (the  surface)  to  r=®.  For  a fixed  energy  E there  is  only  a sincle 

2 

curve  to  analyze,  such  as  that  of  each  example  in  Fig.  2,  which  is  an  r-J 
projection  of  the  phase  space. 

Figure  2 shows  3 examples  illustrating  the  use  of  the  turning-point 
method,  which  depends  on  where  and  how  least  values  ("negative  bumps")  occur 
in  the  g-function  of  Eq.  (59).  The  least  values  may  or  may  not  be  analytic 
minima.  They  correspond  to  "positive  Dumps"  (representing  centrifugal 
barriers)  in  the  effective  potential.  In  these  examples  E-q  is  everywnere 
positive,  but  $ is  otherwise  arbitrary  and  in  particular  may  be  nonmonotonic. 
All  physically  possible  orbits  are  horizontal  lines  constrained  to  remain 
below  the  function  g=r^(E-$).  Orbits  terminating  on  the  g-function  have 
turning-points  there. 

We  consider  Case  (A)  in  detail,  and  summarize  all  3 cases  in  Table  .. 
There  is  one  principal  least  value  (minimum)  in  the  g-function,  at  r=a , and 
no  secondary  least  values.  (This  corresponds  to  an  effective-potential 
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TABLE  i. 

RANGES  OF  r AND  J2  ASSOCIATED  WITH  FIG.  2 

2 ? 2 

J-]  » ^2  ’ ^3  = uPPer  ^ "t  of  range  of  Type-1,  Tyoe-2,  Type-3  orbits 


Case  (A) 

Ranqe  of  r 

(l»a) 

(a,») 

i! 

9a 

9a 

£ 

9 

9 

Case  (B) 

0,a)* 

9b 

- 

9l=9a 

(a.b) 

9b 

- 

9 

(b,® ) 

9b 

g 

- 

2 

(*Type-4  orbits  for  g^<J  <g) 

Case  (C) 

0 ,a) 

9a 

- 

9 

(a.b) 

9a 

9 

- 

(b,c)* 

9a 

9C 

- 

(c,“) 

9a 

9 

- 

2 

(*Type-4  orbits  for  gc<J  <g) 

-32- 


function  with  a single  maximum,  at  r=a.)  At  all  positions  r,  the  range  of 

Type-1  orbits  is  (o,g  ),  where  g Hg(r  )nj,  , and  J is  the  upper  limit  of 

for  Type-1  orbits  in  accord  with  Eq.  (54),  For  r between  r=l  and  r=a, 

2 2 9 2 

Type-3  orbits  occur  for  J in  the  range  J,  so  that  J,  -g;  there  are 

^2^22 

no  Type-2  orbits.  For  r>a,  Type-2  orbits  occur  for  J in  the  range  J,  <o  <g, 

2 1 
so  that  J2“=g;  there  are  no  Type-3  orbits.  The  foregoing  may  be  summarized 

as  in  Table  1 . 

Cases  (B)  and  (C)  nave  secondary  minima  (tnat  of  Case  (B)  is  a least 

value  occurring  at  r=l),  in  addition  to  tne  principal  minimum.  In  all 

cases,  Type-3  orbits  can  occur  only  to  tne  left,  and  Type-2  only  to  tne 

right,  of  tne  principal  minimum.  In  Case  (B),  for  r between  1 and  a,  tne 

upper  limit  for  Type-3  orbits  is  determined  by  the  secondary  minimum  g.  at 

1 2 

r=l  rather  than  by  the  local  value  of  g;  in  this  range  of  r,  values  of  J 
above  g-,  but  less  than  g are  associated  witn  closed  orbits  (Type-4),  as 
noted  by  the  asterisked  note  in  the  table.  Similarly,  in  Case  (C),  for  r 


between  b and  c,  the  upper  limit  for  Type-2  orbits  is  determined  by  tne 

secondary  minimum  g at  r=c  rather  than  by  the  local  value  of  g;  in  this 
c 0 

range  of  r,  values  of  0 above  g but  less  than  g are  associated  witn 

c 

closed  orbits  (Type-4;,  notea  in  iable  1.  A blank  entry  in  tne  J^-c'lumn 

of  Table  1 means  that  there  are  only  Type-1  orbits  for  the  ambient  particles, 

2 

while  a blank  entry  in  the  -column  means  that  there  are  only  Type-1  orbits 
for  the  emitted  particles.  In  either  of  these  cases,  Eq.  (56)  applies 
rather  tnan  Eq.  (55). 

It  is  evident  from  Fig.  2 that  the  existence  of  closed  orbits  is  asso- 
ciated witn  the  existence  of  secondary  minima  in  the  turning-point  function 

(g). 

The  cases  vjhere  g becomes  negative  in  one  or  more  ranges  of  r will  be 

treated  later. 

having  determined  J-  , (Type-2  orbits),  ana  (Tyne-3  orbits)  in 
Table  1 for  the  example  cases  of  Fig.  2,  we  next  show  how  they  are  used  in 
t.  e monoer.ergetic  density  and  flux  factors  of  Secs.  4 and  5. 

Thus,  we  may  now  write  down  tne  monoenergetic  formulas  for  density  ana 
flux,  for  ambient  and  emitted  particles,  for  the  sphere  and  cylinder.  Tnese 
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formulas  are  given  by  Eqs.  (28),  (43),  (35),  and  (47),  for  tne  spherical 
density,  spherical  flux,  cylindrical  density,  and  cylindrical  flux,  respec- 
tively. In  tne  following  equations,  we  set  m=l  if  there  are  only  Type-1 
orbits.  Otherwise,  we  set  m=2  if  the  particles  are  anbient  and  there  are 
Type-2  orbits,  and  we  set  m=3  if  the  particles  are  emitted  and  there  are 
Type-3  orbits.  Moreover,  the  factor  C is  unity  for  ambient  particles  and 
2 for  emitted  particles. 


Density,  Sphere 

Mns(E)=  § j 5-d[E-<H2/r2]  2 


= - 2v£-*-J m2/r2  + /E-o-J^/r2! 


(60) 


Flux,  Sphere 

M-  (E)  = f’Vd(J2/r2)=J  2/r2=J,2 

Jo 

since  the  flux  of  interest  is  at  r=l . 


(61) 


Density,  Cylinder 

H (E)=-  6*d{sin"‘[J2/r2(E-4))] 

in-  w Jq 


4i2sin''/  2/r2(E-<j))  -sin"  !/j,  2/r2(E-q) 
77  m > 


(62) 


Flux,  Cylinder 


ye).  - M,2 

^ I/, 


(63) 


Note  that  the  flux  formulas  are  tne  same  for  emitted  and  ambient  particles. 
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As  a well-known  example  of  the  use  of  the  foregoing  formulas,  we  apply 
them  to  the  problem  considered  by  Bernstein  and  Rabinowitz  (1959),  namely, 
a spherical  prooe  with  monoenergetic  ambient  attracted  particles.  Tne 
potential  falls  off  rapidly  and  monotonical ly,  and  the  turning-point  func- 
tion has  a form  similar  to  that  of  Case  (A)  in  Fig.  2,  namely,  a simple 
minimum.  (This  corresponds  to  a simple  maximi'm  in  the  effective-potential 
function  studied  by  Bernstein  and  Rabir.owitz. ) Using  Taole  1 and  Eqs.  (60)- 
(63),  we  obtain  the  following  results. 

Density,  Sphere  MnS(E)=: 

4 [/E^$  - A-t-g  /r2]  for  r in  (l,a) 

a 

(64) 

2 [/E-$  + /E-<>-g  /r*"]  for  r in  (a,®) 


Flux,  Spnere 


(E)=ga 


(65) 


Density,  Cylinder  Mnc(E)=: 


— sin"  1 v4a/r^(E-<>)  for  r in  (l,a) 

7T  a 

1 -1  / 2 

1 sin  /g./r  (E-p)  for  r in  (a,~) 


(66; 


Flux,  Cylinder 

Mjs(£)=^7  (67) 

The  foregoing  equations  for  density  and  flux,  Eqs.  ( 64 )- ( 67 ) , are  seen  to 
be  identical  except  for  notation  to  these  derived  by  Bernstein  ana  Rabinowitz 
on  the  basis  of  the  effective-potential  formulation.  The  minimum  of  g occurs 
at  r=a,  which  is  the  so-called  "absorption  radius,"  witnin  whicn  only  Type-1 
(one-way)  orbits  contribute  to  the  density  (and  flux). 
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Another  example  which  is  even  simpler  is  the  ideal  Langmuir  limit 
("orbit-limited") , in  which  the  potential  falls  off  less  rapidly  than  the 
inverse-radius-squared.  In  this  case,  the  g-function  for  an  attractive 
potential  rises  monotonically  from  ns  value  £-6  at  the  surface  (not  shown 
in  Fig.  2).  ihen  J-j  =E-s>0  and  J2  =g  for  ambient  particles.  Tne  spherical 
and  cylindrical  fluxes  are  £-yQ  ana  /E-v0,  respectively;  here  we  have  used 
Eqs.  (65)  and  (67),  respectively , with  g replaced  by  E-4>  . This  value  of 
ga  can  be  used  in  the  second  member  of  cq.  (64)  or  Eq.  (66)  to  yield  tne 
orbit-limited  spherical  and  cylinarical  density,  respectively.  For  zero 
potential,  we  have  g =E  and  4>=0,  so  that  we  obtain  the  familiar  for.  _ 

/£(!  + -1  /r^ )/2  and  1-ir  sin”  (1/r)  for  the  spherical  and  cylindrical 

neutral  density,  respectively.  (For  the  purely  monoenergeti c case  E is 
replaced  by  unity;  otherwise,  with  the  factor  /E  the  integral  over  E gives 
unity.) 

Finally,  a numerical  algorithm  suitable  for  the  computer  to  evaluate 
2 2 2 

J-,  , J2  , and  J3  is  the  following.  We  define  g..  at  a set  of  discrete  grid 
points  at  radii  r. , where  the  potentials  are  4^ , by  g^=r..2(E-C:, ) . Here, 
i=l  denotes  the  surface  radius.  Then  the  computer  is  asked  to  find  three 
types  of  least  values  of  the  g^ -array,  namely: 


2 

J-j  = least  {g.j},  all  i 
2 

J2  = least  {g,},  all  j>i 
2 

J^  = least  { g . } , all  j<i 

^ J 


(68) 


ihe  case  where  g is  negative  over  any  range  of  r,  as  would  occur  at  low 

2 

values  of  E with  a repulsive  potential,  is  treated  by  setting  J-  :0.  In 
addition,  for  ambient  particles  we  set  Mr(E)=0  at  r if  g vanishes  anywhere 
between  r and  infinity;  w.nile  for  emitted  particles  we  set  M (E)=0  at  r if 
g vanishes  anywhere  between  r and  the  surface. 


Whereas  the  program  of  Laframboise  (1966)  handles  no  more  than  two 
bumps  in  the  effective  potential,  the  present  program  (PARKSS  - See  Appen- 
dix C)  can  treat  any  number  of  bumps  in  tne  turning-point  function.  Thus, 
the  method  of  this  report  is  completely  general,  and  applies  not  only  to 
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7.  ENERGY  QUADRATURES 

In  order  to  evaluate  numerically  the  integrals  over  energy  in  Secs.  4 
and  5,  one  may  employ  quadrature  formulas  in  the  form: 


r 

J E 


K 

v 


mi  n 


exp(-E)dE*M(E)=  V C.  *M(E  ) 

K=]  K K 


(69) 


where  M(E)  denotes  tne  monoenergetic  factor  dependent  upon  energy  E,  and 
the  lower  limit  En-n  may  De  negative,  in  the  sum  there  are  K coefficients 
C ^ and  abscissas  (energies)  , with  k=l,2,...,K.  The  coefficients  _ .c 

abscissas  depend  on  the  order  K and  nature  of  the  quadrature  scheme  cnosen; 
those  used  in  the  present  work  are  derived  in  Appendix  A. 

This  section  completes  the  analytical  oasis  for  a computer  program. 
Numerical  results  are  presented  next. 
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8.  SAMPLE  RESulTS  AND  DISCUSSION 

A computer  program  named  PARKSS  (for  Parker  Spherical  Sheath)  has 
been  developed  to  implement  the  analysis  of  this  report  pertaining  to  the 
spnerical  problem.  This  section  presents  the  results  of  a number  of  solu- 
tions. The  numerical  parameters  (FORTRAN  input  in  PARKSS  code)  include 
(a)  RB  (=rg)  = outer  grid  radius  di video  by  tne  probe  radius,  (b)  NMP  = 
the  number  of  grid  points  including  the  first  and  last  in  the  interval 
(1,  rg),  and  (c)  ME  = one-half  tne  number  of  energies  in  tne  semi-infinite 
range.  (There  are  also  2 (NMP-1)  additional  values  of  energy  in  the  f 
range  - see  Appendix  A.)  The  3 parameters  RB,  NMP,  and  ME  must  be  large 
enougn  to  ensure  that  the  solution  becomes  insensitive  to  their  values. 

One  also  needs  to  assign  the  value  of  ALPHA,  the  iteration  parameter  con- 
trolling the  mixing  of  successive  iterates. 

For  almost  all  of  the  runs,  RB  was  set  equal  to  2,  NMP  was  21,  ME  was 
5 (resulting  in  50  energy-values ) , and  ALPHA  was  0.1.  Except  for  emission- 
alone  cases  (for  which  the  Dirichlet  condition  $=0  was  used),  tne  inverse- 
square-law  floating  boundary  condition  was  employed  at  tne  outermost  grid 
point  in  all  cases.  The  results  will  be  discussed  with  reference  to  spe- 
cific figures  and  tables. 

A.  Figure  3 and  Tables  2-4 

Figure  3 shows  tnree  sets  of  dimensionless  potential  (<j>)  profiles: 
ambient  plasma  alone  with  no  emission  (curves  labelled  "A"  on  left  side  of 
Fig.  5/,  emitted  electrons  alone  witn  no  ambient  plasma  (curves  labelled 
"E"  on  left  side  of  Fig.  3);  and  both  ambient  and  emitted  (unlabelled 
curves  on  right  side  of  Fig.  3).  Tne  profiles  in  Fig.  3 are  associated 
witn  five  values  of  the  surface  potential  ($,),  namely,  0.8,  0.4,  0.,  -0.4, 
and  -0.8,  corresponding  to  different  vertical  positions  in  the  figure  as 
indicated  by  the  encircled  numbers.  Both  the  ambient  and  emitted  thermal 
fluxes  (Maxwellian)  are  assumed  equal,  corresponding  to  equal  ambient  and 
photoelectron  fluxes,  a condition  which  can  occur  in  the  plasmasphere  (Grard 
et  al,  1973).  In  particular,  the  temperatures  of  the  plasma  ions  and  elec- 
trons (T^,Tg),  of  tne  emitted  electrons  (T^),  are  assumed  to  be  all  equal 
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to  one  electron-volt;  and  tne  density  o plasma  ions  and  electrons  (N  ) is 

equal  to  the  emission  density  of  photoelectrons  (N  1,  and  both  are  taken 
3 s o 

equal  to  400/cm  . Thus,  the  ambient  and  emitted  thermal  currents  are  both 
-9  ? 

about  2x10  amp/cm',  this  being  approximately  the  photoemission  current. 
Thus,  for  a sphere  radius  of  100cm,  tne  Debye  number  is  0.37. 

Referring  to  the  pure-ambient  profiles  (A  in  Fig.  3),  we  see  that  <p  is 
always  monotone  ar.d  the  profile  for  surface  potential  is  the  negative  of 
the  profile  for  Adding  emission  pushes  the  ^-profiles  downward  toward 

negative  values,  and  tne  profiles  are  then  possibly  nonmonotonic.  For  a. 
large  and  positive,  <p  falls  off  more  rapidly  while  for  a large  and  .cive, 
<j>  falls  off  more  slowly,  than  in  the  corresponding  pure-ambient  case. 

The  pure-emission  profiles  (E  in  Fig.  3)  are  similar  to  the  combined- 
case  profiles  (ambient  plus  emitted),  and  DOtn  have  the  following  features: 

For  $ large  ana  positive,  <j>  is  monotone  downward,  essentially  out  to 
infinity.  However,  there  is  always  a negative  minimum  (potential  "barrier"), 
which  moves  inward  from  infinity  and  becomes  deeper  as  b0  becomes  less  posi- 
tive. Numerically,  it  is  difficult  to  "see"  this  barrier  until  <t>  drops 

o 

below  about  unity. 

For  <s  -0,  we  have  a "oure"  barrier,  of  height  (or  depth)  about  3 units 
for  emission  only  and  under  2 units  for  the  combined  case. 

There  is  a critical  negative  value  of  <j>  such  that  the  field  vanishes 
at  the  surface.  This  is  associated  with  the  barrier  having  arrived  at  the 
surface.  The  critical  value  is  about  -0.3  for  emission  only  and  about  -0.4 
for  the  combined  case.  For  <jq  above  the  critical  value,  the  surface  elec- 
trons are  pushed  inward.  Whereas  for  $ , below  tne  critical  value,  the  sur- 
face electrons  are  pushed  outward;  in  this  case,  the  potential  is  monotone 
upward  and  electrons  emitted  from  the  surface  "roll  downhill."  In  a sense, 
one  may  tnink  of  the  barrier  as  a virtuai  one  wnich  has  moved  to  within  the 
surface. 

Comparing  the  ambient-alone  (A)  and  emi tted-alone  (E)  prof. leu  with  tne 
combined-case  prof  iles,  the  latter  appears  roughly  to  be  ap-.roxi;  ;abk-  by 
simply  averaging  the  A and  l nrofiles. 
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Table  2 presents  dimensionless  fluxes  of  plasma  ions  (j . and  j . ) , 
and  plasma  electrons  (jg)  as  functions  of  tne  surface  potential  ($Q),  in 
the  Maxwellian  ambient-only  case.  Here: 

l 

3.J  = ion  flux  normalized  to  ambient-ion  thermal  flux 

= ion  flux  normalized  to  ambient-electron  thermal  flux, 
assuming  hydrogen  ions  and  mass  ratio  M/m=1336 

je  = electron  flux  normalized  to  ambient-electron  thermal  flux 
Note  that  j.'  (<j>  ) = j.  (-<*>),  and  that  equilibrium  (j.=i  ) occurs  for 

1 O c 0 * 0 

$ less  tnan  -2. 

Table  3 presents  dimensionless  escapi ng-electron  flux  (j  ) as  a func- 
tion of  the  surface  potential  ($  ),  in  the  Maxwellian  emission-onl v case, 
i is  normalized  by  the  emitted  flux.  As  $ decreases  from  +2,  j.  increases 

S 0 3 

until  it  becomes  unity  (all  electrons  escaping)  when  passes  the  critical 
(zcro-surface-field)  value,  which  is  at  about  -0.8.  Also  shown  is  0 , 

Mil  Pi 

which  denotes  the  (barrier)  potential  minimum,  while  r denotes  the  posi- 

111 

tion  of  the  minimum  in  sphere  radii.  One  can  see  that  the  minimum  moves 
inward  and  becomes  deeper  as  decreases  from  positive  toward  negative 
values.  At  <>^=0,  tne  minimum  is  -0.28.  There  is  no  equilibrium  surface 
potential  for  Maxwellian  emission-only,  as  there  would  be  for  the  monoener- 
getic  emission-only  case. 

i 

Table  4 presents  dimensionless  fluxes,  of  plasma  ions  (j,  and  ),  of 

plasma  electrons  (j  ),  and  of  escaping  photoelectron  flux  (j  ),  as  functions 

of  the  surface  potential  ($Q),  in  the  combined  Cuse.  Here,  j.,  jfi,  ana  j 

are  all  normalized  by  the  emitted  ."lux  (which  in  this  case  is  the  sa  :e  as 

tne  ambient-electron  thermal  flux),  j has  the  same  definition  as  in  Table 

2.  Table  4 is  qualitatively  similar  to  Table  3.  i increases  as 

s o 

decreases  until  $ passes  the  critical  value,  which  is  at  about  -0.2,  as 

U 

compared  with  -0.8  in  the  emission-only  case.  At  <>  =0,  tne  minimum  is 
-0.12,  as  compared  witn  -0.28  in  the  emission-only  case.  However,  the 
radial  position  of  the  minimum  is  unchanged,  at  1.33.  Tne  equilibrium 
surface  potential  is  at  slightly  above  zero,  compared  with  less  tnan 
-2  in  the  ambient-only  case. 


The  numbers  of  iterations  required  to  achieve  convergence  of  the  Doten- 

5 

tial-iterates  to  one  part  in  10  was  50-60  for  the  ambient-alone,  70-100 
for  the  emitted-alone, and  40-100  for  the  combined  cases. 

B . Figures  4 and  5,  and  Table  5 

Figures  4 and  5 are  for  a combined  problem  with  both  ambient  plasma  anc 
emitted  electrons,  but  where  the  emission  is  monoenergetic  instead  of  Max- 
wellian. The  plasma  temperatures  are  Tj=T  =5eV,  while  the  emission  tempera- 

1 e 3 

ture  is  Ts=leV;  the  plasma  and  emission  densities  are  N0=l/cm  ar.d  7‘1/cm 

respectively.  The  plasma  conditions  are  typical  of  those  in  the  solar  wind 

(Grard  et  al , 1973),  while  the  emitted  current  corresponds  to  photoemission, 
-9  2 

about  10  amp/cm  . This  is  a case  of  "strong  emission,"  defined  by  N /Ts 
»N0/Te,  that  is, the  emission  flux  (j  ) is  much  larger  tnan  the  ambient- 
electron  tnerinal  flux  (jg0).  For  r =100cm,  tne  emission  Debye  number  is 
Xq^O.51,  while  tne  ambient-electron  Deuye  numoer  is  Xqq=16.6. 

In  Table  5,  *t>  is  the  dimensionless  surface  potential  based  on  T , 

whereas  the  dimensionless  fluxes  of  ions  (j.),  of  plasma  electrons  (i  ),  and 

of  escaping  photoelectrons  (j  ),  are  normalized  with  respect  to  tne  emission 

flux.  Qualitatively,  this  table  is  roughly  similar  to  Table  4,  but  there 

are  some  striking  differences.  The  equilibrium  is  very  slightly  below 

unity.  For  tnis  case  of  strong  monoenergetic  emission  it  turns  out  to  oe 

very  difficult  to  obtain  numerical  solutions  very  close  to  ? -1 , namely, 

in  the  range  between  the  entries  <a  =0.99  and  c =1.075  of  the  table. 

o o 

For  <5  above  unity: 

0 

(a)  All  emitted  electrons  return  to  tne  surface,  and  are 
contained  within  a racius  r(n.=Q). 

J 

(b)  There  is  no  minimum;  tne  potential  is  scnctly  monotone. 

(c)  The  position  beyond  which  there  are  no  emitted  electrons, 

denoted  by  r(ns=0),  within  which  all  emitted  electrons  are 

contained,  moves  inward  toward  the  surface  witn  increasing 

<p  above  unity, 
o J 


I 


As  <?0  drops  below  unity,  tnere  is  a minimum  which  moves  inward  from 
infinity  ana  becomes  ceeper.  The  minimum  attains  a aepth  of  -.40  at  <?o-0, 
and  the  condition  of  zero  field  occurs  at  between  $ =-.5  and  $ =-l. 

Figure  4,  where  $o=0.99  (close  to  tne  equilibrium  value),  shows  the 
potential  profile  (dashed  curve  labelled  “<j>")  and  the  profiles  for  dimen- 
sionless plasma- ion  aensity  (n.),  plasma-electron  density  (n  ),  ana  emitted- 
electron  density  (n^).  The  density  profiles  are  the  solid  curves.  The 
emi tted-electron  density  is  1.99  at  tne  surface,  since  almost  all  emitted 
electrons  return.  The  plasma-ion  and  plasma-electron  densities  remain 
close  to  one  another,  cross  over,  and  essentially  approximate  to  the  ne^ur& 
density  (density  at  zero  potential),  shown  as  the  dotted  curve  labelled 
“NEUTRAL".  Hence,  the  ambient  plasma  seems  to  contribute  little  to  the 
space  charge. 

Tne  potential  profile  in  Fig.  4 actually  hua  a minimum  of  -0.0072  at 
r=1.95  (Table  5),  but  this  is  negligible  in  tne  scale  of  the  figure. 

Figure  5 shows  the  same  profiles  as  Fig.  4,  but  for  the  larger  positive 
surface  potential  <j>  =+2.0.  This  is  in  the  regime  where  all  emittea  elec- 
trons return  (ns=2.0  at  the  surface)  ana  are  contained  within  tne  radius 
r=1.42.  This  datum  is  indicated  in  the  second  row  of  entries  in  Table  5 
The  interesting  feature  is  that  n$(r)  approaches  zero  with  infinite  slope, 
at  the  position  where  $=1.0  (where  the  electron  speed  vanishes).  In  this 
case,  n^  and  are  somewhat  more  separated  than  in  tne  $Q=.99  case.  This 
case  was  used  to  make  comparisons  with  a time-simulation  calculation  in 
which  $0=+2  was  found  to  occur  during  a very  slow  transient  following  a 
rapid  or.e  (P.  Rothwell , private  communication,  1976). 

Runs  were  also  made  (not  shown)  in  which  the  emission  was  Maxwellian 
instead  of  nonoenergetic,  with  qualitatively  similar  results.  One  of  the 
cnanges  which  occur,  however,  is  that  finite  but  negligible  potential 
occur  at  large  distances  for  large  positive  $ . 

This  problem  is  related  to  that  of  Schroder  (1973)  in  being  appropriate 
to  a spherically-symmetric  pnotoemissi ve  space  probe  in  the  solar  wind. 
Schr'dder  computed  a number  of  potential  and  density  distributions,  corres- 
ponding to  different  surface  potentials  and  different  plas  a densities,  for 


-43- 


r — i 


i 

a posi ti vely-charged  probe.  The  computational  method  used  is  claimed  to 
follow  that  of  Laframooise  (1966)  but  no  details  are  given.  (The  Lafram- 
boise  method  is  designed  for  monotonic  potentials  only.)  Schroder's 
results  show  the  features  indicated  in  Tables  4 and  5,  ana  Fig.  3,  in  that 
the  potential  minimum  moves  inward  and  oecomes  deeper  as  tne  positive  surface 
potential  decreases.  However,  whereas  Scnroaer's  curves  for  electron  den- 
sity move  Rionotonicdily  downward  wion  increasing  bQ>  it  is  found  in  present 
calculations  that  the  density  at  a given  r rises  to  a maximum  and  then  falls 
off.  Schroder  also  implies  that  the  electron  density  falls  off  wi ir 
inverse-square  cf  the  radius  at  large  radii.  This  is  not  found  to  be  valid 
here  for  positive  surface  potentials,  but  it  is  valid  for  negative  poten- 
tials and  cold  emission,  as  in  the  next  problem. 

The  numbers  of  iterations  required  for  individual  solutions  of  this 

* 

problem  were  in  the  range  50-100. 

C.  Figure  6 

The  problem  of  Fig.  6 again  a combined  ambient-emitted  one,  but  is 
complementary  to  the  previcus  problem  in  that  the  emission  here  is  "wear." 
rather  tnan  strong.  That  is,  the  emission  flux  (j  _)  is  much  smaller  tnan 

SO 

the  ambient-electron  t.nermal  flux  (j  ).  floreover,  the  ambient-electron 
temperature  and  density  are  Tfi- 1 eV  and  NQ- tOO/cm"’,  while  tne  emi ttea-electron 
temperature  T$  is  vanishingly  small  in  toe  sense  that  T,./7e  is  infinitesimal. 
Thus,  the  emission  is  ,!cold"  as  well  as  weak.  The  ambient-electron  Debye 
number  is  0.743  for  a sphere  radius  of  100  cm.  It  is  assumed  for  the  pur- 
poses of  Fig.  6 that  the  emission  flux  is  one-tenth  of  the  ambient-electron 
thermal  flux  (j  /j  =0.1),  and  that  the  ions  are  unperturbed.  Thus,  tne 

S O 60 

equilibrium  potential  for  this  problem  is  *o=-2.3.  (The  ion  motions  were 
included  in  similar  problems  not  shown  here,  but  were  excluded  for  th;r 
problem  to  allow  a comparison  with  an  independent  time-simulation  calcula- 
tion (A.  Wilson,  private  communication,  1976).  Including  the  ion  motif" 
moves  the  equilibrium  potential  up  to  about  -1.8,  but  tne  qualitative 
behavior  of  tne  curves  is  the  same.)  The  assumed  valu<  s of  tt  ■.  oar, meters 
can  occur  in  the  sol«r  wind,  out  with  reduceu  solar  illumination,  say  at 
dawn  or  dusk. 
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In  Fig.  6 the  ^-distribution  at  the  equilibrium  potential  is  shown  as 
a dashed  curve  U is  the  potential  normalized  by  kTe/e).  The  solid  curves 
are  for  both  types  of  electron  density,  with  N /N  denoting  the  dimensionless 
plasma-electron  density,  and  Ns/NQ  demoting  the  emi tted-electror.  density 
normalized  Dy  N . The  emi tted-electron  density  due  to  cold  emission  may  be 
shown  (Appendix  B)  to  be  given  by  the  formula 


Ns(r) 


_J fso  1 

2/7  r2  ^eo  ^'-Oo 


(70) 


which  corresponds  to  electrons  rolling  downhill.  The  numerical  factor 
becomes  1/20/7  for  J so/ J eo=0 • ^ - This  formula  also  indicates  that  Ns/N0 
becomes  infinite  at  the  surface,  which  is  consistent  with  tne  assumption 
of  cold  emission.  Asymptotically,  the  density  crops  as  r”2. 

About  70  iterations  were  required  for  this  problem. 

0 - Figure  7 

Figure  7 shows  the  effects  addin'  surface  emission  to  a probe  problem 
treated  earlier  by  the  author  (Parker,  1975)  which  involved  only  the  ambient 
plasma.  In  the  earlier  problem,  the  dimensionless  potential  and  Debye  number 
were  <j>o=-10  and  X^=0 . 1 , respectively.  In  Fig.  7 the  ambient-only  curves  of 
potential  (9)  and  density  (n^,  nft)  are  labelled  "A",  to  be  compared  with  the 
unlabelled  corresponding  curves  including  emission  (n  ).  The  emission  density 
and  temperature  are  assumed  to  be  the  same  as  the  ambient  values  (equal  emis- 
sion and  ambient  fluxes),  and  all  3 distributions  are  Maxwellian.  This  situa- 
tion can  occur  in  the  plasmasphere  if  the  spacecraft  is  several  meters  in 
radius.  (To  bring  the  equilibrium  potential  cown  to  -10  would  probably  also 
require  the  use  of  an  on-board  electron  gun.)  The  sane  situation  can  also 
occur  in  the  laboratory  with  a small  Langmuir  probe  emitting  thermal  or 
secondary  electrons.  Tne  curves  labelled  "A"  have  been  oresented  in  an 
earlier  report  (Parker,  1975).  The  changes  due  to  emission  are  character- 
ized by  the  unlabelled  curves,  as  follows: 
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(1)  The  ^-distribution  is  shifted  upward  slightly  (in  the  negative 
direction). 

(2)  The  ion  density  is  raisec  about  10  percent  near  the  surface, 
but  tne  effect  on  the  ions  is  insignificant  Deyond  about  1.4 
sphere  radii. 

(3)  The  plasma-electron  density  is  reduced  at  ail  radii. 

(4)  However,  one  quasi  neutrality  point  beyond  which  tne  positive 
and  negative  charges  balance  (n^n.+r  ) to  within  one 
percent  is  still  at  about  the  same  radius  (apout  2 radii)  r. 
in  tne  ambient-alone  case. 

(5)  The  emitted-electron  density  (nc)  drops  off  very  rapidly  at 

first,  and  asymptotically  with  tne  inverse  square  of  the  radius. 

This  behavior  is  associated  with  electrons  rolling  downhill, 

and  is  characteristic  of  the  "cold  emission"  treated  in  tne 

orevious  problem,  where  ,eo  j>>kT  . 

o 1 s 

With  RB=3,  MMP=41 , and  ALPHA=0.01,  about  500  iterations  were  requireo. 

(By  comparison,  about  300  iterations  were  required  when  RB,  NMP,  and  ALPHA 
were  2,  21,  and  0.02,  respectively.) 

£.  Table  6 

In  an  attempt  to  interoret  data  obtained  by  instrumentation  on  the  ATS-6 
spacecraft,  Whipple  (1976)  hypotnesized  the  existence  of  a potential  barrier 
due  to  secondary  electrons  or  photoelectrons.  in  order  to  study  this  effect, 
Whipple  assumed  a spherical ly-symmetri c model  for  the  sheatn  as  in  tne  pre- 
sent report.  However,  his  method  was  approximate  in  that  (a)  simpi if . cations 
were  assumed  for  the  shapes  of  tne  boundaries  in  velocity  space,  i.e.,  Un- 
critical values  of  angular  momentum  (see  Table  1),  and  (b)  the  method  searched 
only  for  the  most  negative  values  of  <?  such  that  simultaneously  d$/dr;... 

The  latter  procedure  is  not  likely  to  yield  pnysically  valid  solutions  since 
tne  boundary  condition  at  infinity  is  ignored;  however,  it  may  be  expected 
to  provide  an  upper  bound  for  the  depth  of  the  barrier  potential.  The  for  ,er 
procedure  (simplifications  in  velocity  space)  should  be  justified  by  exact 
theory  such  as  the  present  one. 
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Table  6 is  concerned  with  the  4 days  (Days  198,  1S9,  204,  ana  273} 
defined  in  Wnipple's  paper.  The  ambient  and  emission  data  for  the  4 days 
are  given  by  the  dimensional  quantities  (N^,  N .T-,  T , T, ) in  the  table. 

G i c S 

The  plasma-ion  and  plasma-electron  temperatures  are  not  equal.  For  each 
day,  the  table  shows: 

(1)  The  relative  strength  of  emission  Current  to  plasma  current. 

(2)  The  surface  potential  $ , negative  or  zero  for  all  4 days. 

(3)  Barrier  potential  (minimum)  $ n , according  to  measurement. 

(4)  Barrier  potential  if  any  4 ^ (W),  according  to  Whipple's 
method. 

(5)  Barrier  potential  if  any  o • (P),  according  to  the  exact 
theory  of  this  report  (computer  program  PARKSS). 

(6)  Barrier  position  if  any  rm(W),  accord:.  . to  Whipple's  method. 

(7)  Barrier  position  if  any  r (P),  according  to  PAAKSb. 

The  results  are  as  follows: 

For  Days  198,  199,  and  204,  the  WlrpDle  method  predicts  a barrier, 
while  PARKSS  predicts  either  no  earner  ^Day  198)  or  a smaller  barrier 
(Days  199  and  2C4).  For  Day  273  both  theories  agree  in  predicting  no 
barrier.  In  all  cases,  the  measurements  indicate  the  existence  of  a bar- 
rier and  tne  deptn  of  this  barrier  is  greater  than  any  of  the  *edictec 
values.  The  discrepancy  between  theory  and  measurement  s striking  T"or 
Day  273,  where  the  emission  is  cold  and  weak,  and  where  tne  e itted  elec- 
trons would  be  expected  to  roll  downhill  without  building  up  a barrier. 

The  potential  distributions  for  tne  4 problems  are  not  shown  here.  Put  are 
typified  oy  several  of  tne  profiles  ’in  rig.  3,  for  zero  and  negative  <?  . 

Tne  numbers  of  iterations  required  for  tne  4 problems  ranqea  from  4,; 
to  60. 
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TABLE  2. 


FLUXES:  AMBIENT  ONLY  (MAXWELLIAN) 


5* « 
u 

ji 

2 

.135 

.0032 

2.77 

1 

.37 

.0086 

1.93 

.8 

.45 

.011 

1.76 

.6 

.55 

.013 

1.58 

.4 

.67 

.016 

1.39 

.2 

.82 

.019 

1.20 

0. 

1.00 

.023 

1.00 

-.2 

1.20 

.028 

.82 

-.4 

1.39 

.032 

.67 

-.6 

1.58 

.037 

.55 

-.8 

1.76 

.041 

.45 

-1 

1.93 

.045 

.37 

_2  ^ 

2.77 

.065 

.135 

(a)  Equilibrium  <j>  below  -2. 


I 
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TABlE  3. 


FLUXES: 

o 

0 

EMITTED  ONLY 
An  = .37 

i 

Js 

(MAXWELLIAN / 

6 • (O 
ymm  rn 

2 

.33 

negligible 

1 

.63 

-.013  (1.85) 

.8 

.70 

1 

o 
. 1 

o 

.6 

.76 

-.082  (1.60) 

.4 

.82 

-.14  (1.50) 

.2 

.88 

1 

ro 

o 

-p> 

LTl 

0. 

.92 

-.28  (1.35) 

-.2 

.96 

-.37  (1.25) 

-.4 

.99 

-.49  (1.20) 

-.6 

1.00 

-.62  (1.10) 

1.00 

monotone 

-.99 

1.00 

monotone 

(a)  Surface  field  nearly  zero. 


TABLE  4. 

FLUXES:  AMBIENT  AND  EMITTED  (MAXWELLIAN) 


= 0.37 

L> 


*0 

ji 

■i 

Js 

Wrm> 

2 

.135 

.0032 

2.61 

. 35 

negligible 

1 

.37 

.0086 

1.80 

.66 

negl igible 

.8 

.45 

.011 

1.64 

.74 

HGQ  J 1 J . J I 

.6 

.55 

.013 

1.47 

.81 

-.009  (1.95) 

.4 

.57 

.016 

1.31 

.88 

-.033  (1.75) 

.2 

.82 

.019 

1.15 

.94 

-.065  (1.55) 

o.u> 

1.00 

.023 

.98 

.98 

-.12  (1.35) 

-.2  (b) 

1.20 

.028 

.82 

1.00 

-.22  (1.10) 

-.4 

1.40 

.033 

.67 

1.00 

monotone 

-.6 

1.60 

.037 

.55 

1.00 

monotone 

-.8 

1.79 

.042 

.45 

1.00 

monotone 

-.99 

1.97 

.046 

.37 

1.00 

monotone 

■2. 

2.86 

.067 

. 1 35 

1 . 00 

monotone 

(a)  Equi  nbrium  e> 

0 

slightly  above 

0. 

(d)  Surface  field 

vanishes  for  <J> 

about  -0.2. 
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A 


D. 


FLUXES:  MAXWELLIAN  AMBIENT  AND  MONOENERGETIC  EMITTED 


<T1 

’ Te’  Ts  = 5 

. 5,  1 eV 

• N N = 1 

u S 0 

, 211/cm3;  rQ  = I 

00  cm) 

C.s 

O 

. (a) 
J1 

, (a) 
Je 

. (a ) 

JS 

‘Vsi  rSr  io; 

r(ns=o) 

A 

J 

. U00u68 

. 0085 

0. 

monotone 

1.28 

2 

.000033 

.0074 

0. 

monotone 

1.42 

1.5 

.000092 

.0069 

n 

u , 

monotone 

1.55 

1.075 

.00010 

.0064 

A 

U . 

monotone 

1.80 

>99(o) 

.00010 

.0063 

.011 

-.0072  (1.95) 

- 

.5 

.00011 

.0058 

.58 

-.2.  (1.70) 

- 

0. 

.00012 

.0053 

1.00 

-.40  (1.55) 

- 

_ r 

.00014 

.0048 

i 

-.e9  v . 28 ) 

- 

-1 

.00015 

.0043 

1.00 

monotone 

- 

(a)  Ion 

mass  = 1836 

electron 

masses.  All 

fluxes  normalize 

0 oy  emission 

fl  ex. 

(b)  Equilibrium  9.  slightly  less  tnan  +1.0. 


3 I - 


TABLE  6. 


ATS  6 DATA  AND  THEORETICAL  COMPARISONS 

ambient  density  (cm  ; 

photoelectron  density  vcri  ) 

ambient-ion  temperature  (eV' 

ambient-electron  temperature  (eV) 

photoelectron  temperature  (eV) 

satellite  potential  (volts) 

barrier  potential  (volts) 

barrier  racius  divided  by  ro=5m 

barrier  potential,  Whipple  approximation  (1976) 

barrier  potential,  Parxer  eAcet  theory  (PARK5S) 

oamer  radius,  Whipple  approximation  (1976) 
barrier  radius,  Parxer  exact  theory  (PARKSS) 


Quanti ty 

Day  198 

Day  199 

Day  204 

Dag  273 

N 

0 

0.2 

1.2 

90 

17 

Nso 

10.3 

90 

200 

o.< 

T . 

T 

7 

10 

3.5 

650 

i 

e 

DO 

32 

14.5 

3u0 

s 

4.9 

2 

1 .9 

6 

Relative  strength 
of  pnotoemission* 

strong 

strong 

conpcti ti ve 

W6aK 

-20 

0 

9 

- L 

-2000 

«>..  (measured) 

i.iin 

-60 

- iu 

-10 

-2050 

W“> 

-24.5 

-6.2 

-3.2 

none 

W?) 

none 

-2  2 

-2.5 

none 

r_(W) 

6 

4 

1.5 

none 

rm(P) 

none 

1.5 

1.2 

none 

N /7~  versus  N /T- 
so  s o e 


min 


rm 

‘min 
rm(W) 

li] 

r (P) 


J , = 

min’ ' 

= 


APPENDIX  A. 


ENERGY  QUADRATURE  FORMULAS 

The  energy  quadrature  formulas  used  in  this  work  have  the  form,  appro- 
priate  to  Maxwellian  energy  attributions: 


[const) 


e/?(-E)dE*M(E)  - 


'min 


<=1 


'k  ‘'ivV 


( A— 1 ) 


where  "const"  denotes  a constant  multiplier,  and  where  M(E)  is  tne  mono- 
energetic  factor  depending  on  the  energy  E resulting  from  integration  over 

9 

o . Ck  and  Ek  denote  the  coefficients  ana  abscissas  respectively  of  the 
K-th-order  quadrature  formula.  The  coefficients  are  evaluated  by  PARKbS 
as  follows: 


Given  a potential  distribution  defined  at  a set  of  radial  grid  points 
r, . there  is  an  associated  set  of  discrete  potentials  ^^(t- ).  Arrange 
tne  set  of  (^-values  in  order  of  increasing  value.  Then,  if  E.^  of  Eq. 

(A-l)  is  identified  with  the  least  value-  of  the  set  of  o -values,  and  if  E 

i rax 

is  identified  with  the  greatest  value  of  this  set,  we  may  split  the  integral 
ana  sum  of  Eq.  (A-l)  into  two  parts,  one  corresponding  to  the  finite  range 
(E.,  . Em  ),  and  tne  ether  to  the  semi -infinite  range  (E„,  . <»): 

nil  N hid  A lilaX 


exp(-E)de*M(E) 


'max 


min 


r°° 

exp(-E)d£*M(E) 

max 


(A-2) 


where  the  const  factor  of  Eq.  (A-l)  has  been  omitted. 

First  we  consider  the  seni-finite  range  (seccr.a  integral),  and  then 
tne  finite  range  (first  integral). 


A . 1 Semi -Infinite  Range 

For  the  semi -infinite  range  we  may  employ  quad  riture  formulas  develc,..,-.: 
by  Steen  et  al  (1060),  or  by  Eaframboiso  and  Stauffer  ; lhiV<).  These  ar.- 
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designed  for  the  Maxwellian  case  where  tne  integrand  contains  a Gaussian 
function  as  a weighting  function.  Let  a.,,  H.  denote  an  abscissa-coeffi- 
cient  pair  from  the  data  of  Steen  et  al  (1969).  Then  one  transforms  the 
semi -infinite  range  integral  to  ^ 


k 


max 


exp(-E)d£*M(E) 


= exp(-E  ) exp(-u)dU'M(li+E  ) 

liidx  , max 

0 


N 

u 

k=l 


2Hkak‘B(ak 


"max' 


(A-3) 


Here  we  have  formed  the  k-th  coefficient  as  C, =2n.  a,  , and  the  associated 

? K K K 

,k-tn  energy  abscissa  as  E,„=ai~  +E_, v.  The  coefficients  H,  and  abscissas  a,. 

K K.  litaX  K K 

may  also  be  obtained  from  the  "One-Dimensional"  table  of  Lafranboise  and 
Stauffer  if  one  multiplies  their  coefficients  by  /if  and  their  abscissas  by 
unity.  (There  is  also  the  option  of  using  the  Laframboi se-Stauffer  "Two- 
Dimensional"  or  "Three-Dimensional"  abscissas  and  coefficients  as  they  sug- 
gest for  a cylindrical  or  spherical  problem,  respectively.  However,  it  is 
the  author's  opinion  based  cn  experience  with  ootr,  methods  that  this  yields 
no  significant  gain  in  practice.) 


A. 2 Finite  Range 

The  finite  range  of  energies  (E..,  , E„v)  defined  by  the  set  cf  ordered 

nun  iTicix 

values  of  4* . consists  of  a number  of  unequal  energy  intervals.  The  number 
of  sjcii  intervals  is  equal  to  the  number  of  grid  points  minus  one.  Consider 
one  of  those  intervals,  and  assume  its  energy  range  is  (A, 8).  Then  its 
ccr.triDution  to  the  energy  integral  may  be  written 
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-B 


exp(-E)dE>M(t) 

jA 


' - A.  ry 

exp^-A-u^ ) •2udu'M(u‘‘+A) 

Q 


= y Ck'M(Ek=uk“  +A)  (A-4) 

k 

To  obtain  a given  desired  order  of  tne  quadrature  formula,  one  may  obtain 
from  £q.  (A-4)  a sequence  of  moment  equations.  For  Order  2 as  used  in 
PARKSS,  we  have: 


,2\  m 


T I / i\  Cm  \ III  r\  . r*  m , r m 

Im  = exp(-A-u  )*u  •dudu  - C^u-. 

1 r\ 


(A-5) 


where  m-0,1,2,  and  3.  Equations  (A-5)  const! tutu  4 equations  for  4 unknowns: 
u.j , u9,  C-j , and  C2.  Tne  solution  is  as  follows: 


u.j  = b 


-Jfc 


| u^  = b + A?-a 


(A-6) 


[ a = 0.5(I0I3-I,I2)/(I0I2-I 


2 V 
1 } 


a (I2  -I-j I3)/ ( I-J  _1oV 


( A- 7 ) 


j Cj  2a-j exp ( -A ; v I-j “iQu2y 


i v ^ ^ o / 


^2  ~ 2u2@xp(“A ) v 1 -j - ^q^*|  )/ ^2*^1 ) 


(A-8) 


.oere 


-A3- 


I0  = \ erf(^) 
l]  = 0.5(l-exp(A-B)) 

= 0.5( Iq-/B-A  exp(A-B)) 
I3  = I^-0.5(S-A)  exp(A-S) 


(A-9) 


It  may  be  noted  tnat  replacing  B by  infinity  and  A by  zero  in  the  above 
equations  yields  the  abscissas  and  coefficients  of  Steen  et  al  for  Cm  ~ 2. 

For  small  (B-A),  one  may  expand  the  foregoing  equations  to  obtain  the 
lowest  order: 


C-,  = = C = G.5/3-A 

u1  = (1-1//3)C  ( A- 10) 

u2  s (1+1  //3)C 


After  computing  the  u's  and  C 1 s as  above,  one  forms  the  coefficients 
ana  energies  used  in  PARKSS  by: 

E,  = u/  + A 

i i 

- - 2 + £ 
c0  - U0  + A 

(A-ll ) 

C-j  = (const)C-. 

C3  = (const)C2 


wnere  tne  const  factor  is  that  or,  tne  left-hand  side  of  Eq.  (A-l). 

One  nay  choose  an  order  M higher  than  2 for  any  energy  interval,  in 
which  case  the  4 equations,  £q.  (A-5),  are  replaced  by  2M  simultaneous 
equations.  However,  the  additional  computer  tine  requires  may  not  be 
justihioo. 
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APPENDIX  B. 

LIMIT  OF  COLD  EMISSION  FROM  NEGATIVE  PROBE 


We  first  derive  the  expression  for  the  density  of  emitted  particles 
when  the  potential  is  repulsive  U>0)  and  nonotomc.  (Here,  <->  denotes  she 
potential  divided  by  - • kT  /e j ; thus,  this  applies  to  electrons  emitted  from 
a negative  probe  having  dimensionless  potential  with  reference  to  elec- 

U 

trons.)  If  the  emission  is  weak  toe  floating  potential  referred  to  ions 
will  normally  be  negative  (positive  when  referred  to  electrons)  and  t.e 
profile  will  be  monotonic. 

2 

Consider  the  turning-point  function  g r ( E— ; . In  the  above  situation, 

9 is  a oositive  and  monotomcal ly  decreasing  function  of  r,  and  E>o,  so  that 

g also  rises  monotonically.  For  E<o  , we  hav«_  g<0  at  r=l ; hence  there 

0 

can  be  no  emitted  contributions.  For  £>g  we  have  the  least  value  of  q, 

g=E-9  , occurring  at  r=l.  Thus,  there  are  only  Type-1  orbits  (Sec.  6) 

u 2 

occu'iable  by  emitted  electrons,  and  so  we  may  use  Eg.  (5c)  with  J,  =E-p  ■ 

• G 

n view  of  Eqs.  (25)  and  (33)  for  tne  spherical  and  cylindrical  density, 
respectively,  Eqs.  (53)- (56)  yield: 


Mns(E)  = y'T-~  - A -9-(E-9oVr 
(sphere) 


M.  ,(E)  = 

.It 

(cylinder) 


2 

71 


r (E-$) 


(B-l 


Substituting  Eq.  (B-l)  into  Eq.  (25)  (whore  9S  is  replaced  by  and  inte- 
grating from  to  infinity  yields  for  .r.e  spherical  density  due  to  j '‘ax- 
well  iar.  distribution  of  emitted  particles: 


ER(;  ■ 


'?) 


/ o ' r ( 9 . - 
/ 1-1//  ER  / — f/SL 

r - 1 

I 

L 


(i 
\ L 


wr.ere  lR(x)  exp(  x)  *erfox. 


1 


Sir’ll, ar  substitution  ct  Eq.  (3-2)  into  Eq.  (32)  for  tne  cylindrical 
prcblem  results  in: 


so  jo 


~ \ / L 

exp(-u)du‘sin  /u/r‘  (u*-.;,  -i/i) 


(3-4' 


which  cannot  be  expressed  in  terms  of  elementary  functions. 


For  monoenergetic  emission,  we  substitute  1+$  for  E in  Eqs.  (3-1) 


and  (8-2)  so  that 


so 


f~ L ? 

- /]+$  -<p  - -<p  - 1/r 


(3-s ) 


(monoenergetic,  sphere) 


and 


v .2  . -i  / r 

N ' 77  S1n  / 2,,  , 

so  r ( !+<<>-$) 


(8-6) 


ro 


(monoenergetic,  cylinder) 


Now  proceeding  to  the  limit  of  cold  emission,  we  let  9 -$  become  lar  e 


compared  with  unity.  Then  Eq.  ( B— 3 ) is  replaced  by 


(sphere,  Maxwellian) 


so 


(3-7 ) 


Equation  (3-4)  is  replaced  by 


so 


r/q" 

0 ' 


(cylinder,  Maxwellian) 


(8-3) 
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Equation  (B-5)  is  replaced  by 


s 

IT* 

so 


1 


(sphere . monoenerqutic) 


(13-9 1 


And  Eq.  (8-6)  is  replaced  by 


Trr/<{>  -<!> 


(cylinder,  monoenergeticj 


These  limiting  formulas  are  identical  to  tncse  which  would  result  from  the 
assumption  of  purely  radial  motion,  with  no  initial  kinetic  energy.  That 
is,  the  electrons  roll  downhill  after  essentially  starting  from  "rest''. 


APPENDIX  C. 


COMPUTER  PROGRAM  PARKSS 


The  following  FORTRAN  listing  is  for  the  PARK SS  computer  program. 
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